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ABSTRACT 

We present a study of the mass-metallicity (M-Z) relation and H II region physical conditions in a sample of 
20 star-forming galaxies at LO < z < L5 drawn from the DEEP2 Galaxy Redshift Survey. We find a correlation 
between stellar mass and gas-phase oxygen abundance in the sample and compare it with the one observed 
among UV-selected z ^ 2 star-forming galaxies and local objects from the Sloan Digital Sky Survey (SDSS). 
This comparison, based on the same empirical abundance indicator, demonstrates that the zero point of the 
M-Z relationship evolves with redshift, in the sense that galaxies at fixed stellar mass become more metal-rich 
at lower redshift. Measurements of [O III]/H/? and [N II]/Ha emission-line ratios show that, on average, objects 
in the DEEP2 LO < z < L5 sample are significantly offset from the excitation sequence observed in nearby 
H II regions and SDSS emission-line galaxies. In order to fully understand the causes of this offset, which 
is also observed in z 2 star-forming galaxies, we examine in detail the small fraction of SDSS galaxies that 
have similar diagnostic ratios to those of the DEEP2 sample. Some of these galaxies indicate evidence for AGN 
and/or shock activity, which may give rise to their unusual line ratios and contribute to Balmer emission lines 
at the level of 20%. Others indicate no evidence for AGN or shock excitation yet are characterized by higher 
electron densities and temperatures, and therefore interstellar gas pressures, than typical SDSS star-forming 
galaxies of similar stellar mass. These anomalous objects also have higher concentrations and star formation 
rate surface densities, which are directly connected to higher interstellar pressure. Higher star formation rate 
surface densities, interstellar pressures, and H II region ionization parameters may also be common at high 
redshift. These effects must be taken into account when using strong-line indicators to understand the evolution 
of heavy elements in galaxies. When such effects are included, the inferred evolution of the M-Z relation out 
to z ~ 2 is more significant than previous estimates. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

The abundance of heavy elements in the H II regions of 
galaxies reflects the past history of star formation and the ef- 
fects of inflows and outflows of gas. A characterization of 
the evolution of chemical abundances for galaxies of different 
masses is therefore essential to a complete model of galaxy 
formation that includes the physics of baryons (De Lucia et al. 
2004; Finlator & Dave 2007). Important observational con- 
straints for such models come from determining the scaling 
relations at different redshifts among galaxy luminosity, stel- 
lar mass, and metallicity, which, for star-forming galaxies, 
typically consists of the oxygen abundance. However, one 
of the key challenges is to take the observationally measured 
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quantities, i.e. strong, rest-frame optical emission-line ratios, 
and connect them with the physical quantity of interest, i.e. 
oxygen abundance. 

In the local universe, Tremonti et al. (2004) have used a 
sample of ~53,000 emission-line galaxies from the Sloan 
Digital Sky Survey (SDSS) to investigate the luminosity- 
metallicity (L-Z) and mass-metallicity (M-Z) relationships. 
For this sample, metallicities were estimated on the observed 
spectra of several strong emission lines, including [O II] 
AA3726, 3729, H/3, [O III]AA5007, 4959, Ha, [N II] AA6548, 
6584, and [S II] AA6717, 6731. At increasing redshifts, as the 
strong rest-frame optical emission lines shift into the near-IR, 
metallicities are typically based on smaller subsets of strong 
emission lines through the use of empirically calibrated abun- 
dance indicators (e.g. Pettini & Pagel 2004; Pagel et al. 1979). 
Much progress has been made recently in assembling large 
samples of star-forming galaxies with abundance measure- 
ments at both intermediate redshift (Kewley & Dopita 2002; 
SavagUo et al. 2005) and at z > 2 (Erb et al. 2006a). However, 
we have only begun to gather chemical abundance measure- 
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TABLE 1 

Galaxies Observed with Keck II NIRSPEC 



DEEP ID R.A. (J2000) Dec. (J2000) zhc B R I Mb U-B 



42044579 02 30 43.46 00 42 43.60 1.0180 23.22 22.97 22.40 -21.20 0.54 

22046630 16 50 13.83 35 02 01.78 1.0225 23.64 23.02 22.31 -21.37 0.69 

22046748 16 50 14.55 35 02 04.31 1.0241 24.43 23.76 22.90 -20.86 0.86 

42044575 02 30 44.85 00 42 51.33 1.0490 23.08 22.94 22.56 -21.06 0.34 

42010638 02 29 08.74 00 23 28.40 1.3877 22.93 22.85 22.54 -22.12 0.49 

42010637 02 29 08.74 00 23 32.87 1.3882 24.20 23.98 23.72 -20.87 0.44 

42021412 02 3044.55 00 30 50.73 1.3962 24.07 23.74 23.12 -21.91 0.78 

42021652 02 30 44.70 00 30 46.19 1.3984 22.97 22.24 21.32 -24.01 1.01 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes and arcseconds. 



ments for galaxies at z ~ 1-2 (Maier et al. 2006; Shapley 
et al. 2005, hereafter Paper I,). In this work, we continue our 
efforts to fill in the gap of chemical abundance measurements 
during this important redshift range, which hosts the emer- 
gence of the Hubble sequence of disk and elliptical galaxies 
(Dickinson 2000), and the buildup of a significant fraction of 
the stellar mass in the universe (Drory et al. 2005; Dickinson 
et al. 2003) prior to the decline in global star formation rate 
(SFR) density (Madau et al. 1996). 

Chemical abundances for high-redshift galaxies are com- 
monly estimated using locally caUbrated empirical indicators. 
Yet it is crucial to recognize the fact that a considerable frac- 
tion of the z 1 and 2 galaxies with measurements of multiple 
rest-frame optical emission lines do not follow the local exci- 
tation sequence described by nearby H II regions and star- 
forming galaxies in the diagnostic diagram featuring the [O 
III] A5007/H/3 and [N II] A6584/Ha emission-Une ratios (Pa- 
per I; Erb et al. 2006a). On average, the distant galaxies lie 
offset towards higher [O III] A5007/H/3 and [N II] A6584/Ha, 
relative to local galaxies. As discussed in Paper I and Groves 
et al. (2006), several causes may account for this offset, in 
terms of the prevaiUng physical conditions in the H II regions 
of high-redshift galaxies. The relevant conditions are H II re- 
gion electron density, hardness of the ionizing spectrum, ion- 
ization parameter, the effects of shock excitation, and contri- 
butions from an active galactic nucleus (AGN). It is still un- 
clear which of these are most important for determining the 
emergent spectra of high- redshift galaxies. Understanding 
this offset in emission-line ratios is important, not only be- 
cause it provides evidence that physical conditions in the high 
redshift universe are different from the local ones, but also be- 
cause the application of an empirically calibrated abundance 
indicator to a set of H n regions or star-forming galaxies rests 
on the assumption that these objects are similar, on average, 
to those on which the calibration is based. 

In this sense, the current work is also motivated by the in- 
terpretation of the offset in emission-line ratios among distant 
galaxies, and an assessment of the reliability of using local 
abundance calibrations for high-redshift star-forming galax- 
ies. Instead of focusing on high-redshift objects, another ap- 
proach is to study the properties in a class of nearby objects, 
which exhibit similar offset behavior on the emission-line di- 
agnostic diagram. Unravelling the relations between the phys- 
ical conditions and unusual diagnostic line ratios for such ob- 
jects aids the understanding of high-redshift galaxies. The 
SDSS, with its rich set of photometric and spectroscopic in- 
formation, provides an ideal local comparison sample. 

In this paper we expand on the analysis presented in Paper 
I, with an enlarged sample of DEEP2 star-forming galaxies 
observed with NIRSPEC on the Keck II telescope. The larger 



number of DEEP2 objects with near-IR observations enforces 
the conclusions drawn in the previous work. Furthermore, our 
detailed study of nearby SDSS objects with similar emission- 
line diagnostic ratios leads to a clearer physical explanation of 
the observed properties of the DEEP2 galaxies. The DEEP2 
sample, near-IR spectroscopic observations, data reduction, 
and measurements are described in §2. We present the oxygen 
abundances derived from measurements of [O III], H/3, Ha, 
and [N II] emission lines in both individual as well as com- 
posite spectra in §3. The mass-metallicity relationship and 
its evolution through cosmic time are discussed in §4. In §5 
we investigate differences in z ~ 1.0- 1.5 H II region phys- 
ical conditions with respect to local samples by examining 
nearby SDSS galaxies with similar emission-line diagnostic 
ratios. Finally, we summarize our main conclusions in §6. A 
cosmology with f2„ = 0.3, Qa = 0.7, and h = 0.7 is assumed 
throughout. 

2. DEEP2 SAMPLE, OBSERVATIONS, AND DATA REDUCTION 
2.1. DEEP2 Target Sample and Near-IR Spectroscopy 

The high-redshift galaxies presented in this paper are drawn 
from the DEEP2 Galaxy Redshift Survey (hereafter DEEP2; 
Davis et al. 2003; Faber et al. 2007), which contains > 30, 000 
galaxies with high-confidence redshifts at 0.7 < z < 1.5 down 
to a limiting magnitude of Rab = 24.1. The motivation for our 
follow-up near-infrared spectroscopic program, along with 
detailed descriptions of the sample selection, optical and near- 
IR photometry, and spectroscopy are presented in Paper I. 
Only a brief overview is given here. 

The new sample contains four galaxies at z ^ 1.0 and four 
galaxies at z ~ 1.4, which, in combination with the pilot pro- 
gram presented in Paper I, leads to a sample of 20 galaxies 
in total. These galaxies are located in fields 2, 3, and 4 of 
the DEEP2 survey, at 16, 23, and 2 hr right ascension, re- 
spectively. To probe chemical abundances and H II region 
physical conditions, observations of several strong H II region 
emission Unes are required, ideally at least [O II], H/?, [O III], 
Ha, and [N II]. At z > 0.85, however, the only strong H II 
region emission feature contained in the DEEP2 DEIMOS 
spectroscopic data is the [O II] doublet. Therefore, near-IR 
spectroscopy is needed to measure longer wavelength H II 
region emission lines at z > 1. We target two narrow red- 
shift windows within the larger DEEP2 redshift distribution: 
0.96 < z < 1.05 and 1.36 < z < 1.50, within which it is pos- 
sible to measure the full set of H/3, [O III], Ha, and [N II], in 
spite of the bright sky lines and strong atmospheric absorption 
in the near-IR. 

The absolute B magnitude, Mb, and stellar mass estimates 
are given in Tables 1 and 2 for the new objects and plotted as 
red squares in the lower panels of Figure I, together with the 
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Fig. 1 . — DEEP2 color-magnitude and magnitude-stellar mass diagrams. The rest-frame (U—B) vs. Mg color-magnitude diagrams (top) and the Mb vs. stellar 
mass diagrams (bottom) are for DEEP2 galaxies at 0.96 <z< 1-05 (left) and 1.36 < z < 1.50 (right). In each plot, DEEP2 galaxies from both the pilot sample 
presented in Paper I as well as objects with new NIRSPEC observations are indicated with red squares. As shown here, all NIRSPEC targets were drawn from 
the "blue cloud" of the color bimodality. . 



data for the pilot sample. For all the objects, we use the Mb 
values from Willmer et al. (2006) based on optical data alone 
and confirm their good agreement with estimates based on fits 
to the BRIKs SEDs that span through rest-frame J (I) band for 
the objects at z ~ 1 .0 (z ^ 1 .4). Stellar masses for the objects 
in our sample are derived with Ks-hand photometry, follow- 
ing the procedure outlined in Bundy et al. (2005), which as- 
sumes a Chabrier (2003) stellar initial mass function (IMF). 
As discussed in detail in Paper I, the Bundy et al. (2005) stel- 
lar mass modelling technique agrees well with that used by 
Kauffmann et al. (2003b) for SDSS galaxies, based on both 
spectral features and broadband photometry. Stellar masses 
for galaxies in the pilot sample have been updated to reflect 
both the most current DEEP2 near-IR photometric catalog, 
and population synthesis models limiting the stellar popula- 
tion age to be younger than the age of the Universe at z ^ 1 . 
Thus, in a few cases, the stellar masses differ slightly from 
those presented in Table (2) of Paper I. As shown in the lower 
panels of Figure (1), the z ~ 1.4 galaxies in our sample span 
the full range of absolute B luminosities in the DEEP2 survey, 
from Mb ^ -20 to -23, while the smaller set of z 1.0 galax- 
ies happens to cover the faint end of the luminosity function. 
Galaxies in both redshift intervals cover more than an order 
of magnitude in stellar mass and therefore should be able to 
probe an interesting dynamic range. Since our goal was to 
study the emission-line properties of galaxies, all objects in 
our sample lie in the blue component of the observed (U-B) 
color bimodality in the DEEP2 survey, as shown in the upper 
panels of Figure 1 . 

The near-IR spectra were obtained on 2005 September 17 
and 18 with the NIRSPEC spectrograph (McLean et al. 1998) 
on the Keck II telescope. Over the range of redshifts of the 
galaxies presented here, two filter setups are required to mea- 



sure the full set of H/3, [O III], Ha, and [N II]. For objects at 
z 1.4, the NIRSPEC 5 filter (similar to H band) is used to 
observe Ha and [N II], whereas the NIRSPEC 3 filter (similar 
to J band) is used for H/3 and [O III] . For objects at z ^ 1.0, the 
NIRSPEC 3 filter is used to observe Ha and [N II], whereas 
the NIRSPEC 1 filter (AA = 0.95-1.10 fim) is used for H/3 
and [O III]. All targets were observed for 3x900 s in each 
filter with a 0.76 x42 long slit. The spectral resolution de- 
termined from sky lines is 10 for all four NIRSPEC filters 
used here. Photometric conditions and seeing were variable 
throughout both nights, with seeing ranging from 0.5 to 0.7 
in the near-IR. In order to enhance the long-slit observing ef- 
ficiency, we targeted two galaxies simultaneously by placing 
them both on the slit. 

We observed a total of 10 DEEP2 galaxies, successfully 
measuring the fuU set of H/3, [O III], Ha, and [N II] for eight 
out of 10. For the remaining pair, we only detected Ha in 
the H band, but no H/3 nor [O III] in the /-band exposures, in 
which the background in between sky lines was characterized 
by a significantly higher level of continuum than usual. This 
anomalous background is likely due to an increased contri- 
bution from clouds, which may have affected both H- and J- 
band observations of the pair Since a clear measurement was 
not obtained for these two objects, due to variable weather 
conditions, we exclude them from our study. The object, 
42021652, has a double morphology, with one component 
dominated by emission lines with a weak continuum, and an- 
other component dominated by strong continuum, with only 
weak emission lines at roughly the same redshift. The separa- 
tion between the two components on the sky is ^^0.9 , which 
corresponds to ~8 kpc at z = 1.3984, perhaps indicative of a 
merger event. This interpretation is supported by the small ve- 
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TABLE 2 

Emission Lines and Physical Quantities. 



12 + log(0/H) 



DEEP ID 






F[OIII]A5007'' 


Fhc" 


F[NII]A6584'' 


N2'' 


03N2'^ 




SFRHa = 


log(M*/M0)f 


42044579 


1.0180 


2.4±0.4 


5.4±0.3 


12.1±0.3 


2.8±0.3 


8.54±0.18 


8.41±0.14 


0.7 


3 


10.26±0.15 


22046630 


1.0225 


3.7±0.5 


8.7±0.4 


15.2±0.3 


1.8±0.3 


8.37±0.18 


8.32±0.14 


0.8 


3 


10.29±0.16 


22046748 


1.0241 


2.2±0.4 


6.4±0.3 


9.6±0.2 


2.3±0.2 


8.55±0.18 


8.39±0.14 


0.5 


2 


10.28±0.10 


42044575 


1.0490 


7.7±0.7 


22.9±0.3 


23.6±0.3 


3.6±0.3 


8.43±0.18 


8.32±0.14 


1.4 


6 


9.74±0.06 


42010638 


1.3877 


3.8±0.9 


17.3±0.4 


16.6±0.5 


2.8±0.5 


8.46±0.19 


8.27±0.15 


2.0 


9 


10.21±0.06 


42010637 


1.3882 


3.3d=0.9 


3.7±0.3 


6.4±0.4 


2.9±0.4 


8.70±0.18 


8.60±0.15 


0.8 


4 


9.96±0.12 


42021412 


1.3962 


7.1±0.8 


> 2.0 


12.7±0.3 


2.9±0.3 


8.54±0.18 


< 8.70 


1.5 


7 


10.89±0.13 


420216528 


1.3984 


4.6±0.5 


9.7±0.5 


9.5±0.3 


1.2±0.3 


8.38±0.19 


8.33±0.15 


1.1 


5 





'^Emission-line flux and random error in units of 10 ergs s cm . 
''Oxygen abundance deduced from the N2 relationship presented in Pettini & Pagel (2004). 
'^Oxygen abundance deduced from the 03N2 relationship presented in Pettini & Pagel (2004). 
^Hc« limiinosity in units of lO"^ ergs s"'. 

^Star formation rate in units of Mq yr~', calculated from Lho using the calibration of Kennicutt (1998), and divided by a factor of 1.8 to convert to a Chabrier 
(2003) IMF from the Salpeter IMF assumed by Kennicutt (1998). Note that SFRs have not been corrected for dust extinction or aperture effects, which may amoimt 
to a factor of 2 difference (Erb et al. 2006b). 

^Stellar mass and imcertainty estimated using the methods described in Bundy et al. (2005), and assuming a Chabrier (2003) IMF. 

»This object has a double morphology. The separation between the two components is about 0.9 , which corresponds to ~8 kpc at z = 1.3984. We measured 
line fluxes for the emission-line-dominated component, but do not have a robust estimate of the corresponding stellar mass. In the DEEP2 photometry the two 
components were counted as one source, and the resulting stellar-mass estimate has contribution from both components. We therefore do not include this stellar 
mass estimate in otir sample. 



locity difference of Av ~ 1 25 km s"' between the two compo- 
nents. We measure line fluxes for the component dominated 
by emission lines, since it provides a more robust estimate of 
line ratios. Deblended optical and near-infrared magnitudes 
would be required to obtain robust stellar masses for the indi- 
vidual components. However, in the DEEP2 photometry the 
two components were counted as one source since they are 
too close to be deblended and the stellar mass has contribu- 
tion from both of them. For now, we only include flux mea- 
surements of the emission-line component for the diagnostic- 
line-ratio analysis but do not include this object in the mass- 
metallicity studies. A summary of the observations including 
target coordinates, redshifts, and optical and near-IR photom- 
etry is given in Table 1 . 

2.2. Data Reduction and Optimal Background Subtraction 

Data reduction was performed with a similar procedure to 
the one described in Paper I and Erb et al. (2003), with the 
exception of an improved background subtraction method ap- 
plied to the two-dimensional galaxy spectral images (Kel- 
son 2003; Becker 2006, private communication). In the cus- 
tom NIRSPEC long-slit reduction package written by D. G. 
Becker (2006, private communication), optimal background 
subtraction is performed on the unrectified science frames. 
First, a transformation is calculated between CCD (x, y) co- 
ordinates and those of slit position and wavelength, using 
the wavelength-dependent traces of bright standard stars and 
the spatially dependent curves of bright sky lines. Then a 
two-dimensional model of the sky background is constructed 
as a function of slit position and wavelength, using a low- 
order polynomial in the sUt-position (dimension, and a fo-spUne 
function in the wavelength dimension. This two-dimensional 
model is iteratively fit in the differenced frame of adjacent sci- 
ence exposures and subtracted from the unrectified data. Af- 
ter background subtraction, cosmic rays were removed from 
each exposure, which was then rotated, cut out along the slit, 
and rectified. Finally, all background-subtracted, rectified ex- 
posures of a given science target were combined in two di- 
mensions. This new approach to reducing NIRSPEC spec- 



tra results in fewer artifacts around bright sky lines and cos- 
mic rays, which are commonly introduced when rectification 
is performed before sky subtraction and cosmic-ray zapping. 
One-dimensional spectra, along with error spectra, were then 
extracted and flux- calibrated using A-star observations, ac- 
cording to the procedure described in Paper I and Erb et al. 
(2003). 

2.3. Measurements and Physical Quantities 

One-dimensional, flux-calibrated NIRSPEC spectra along 
with the 1 fT error spectra of galaxies in our new sample are 
shown in Figures 2 and 3. Emission-line fluxes and uncertain- 
ties measured from the one-dimensional spectra are given in 
Table 2. Ha and [N 11] A6584 emission-line fluxes were de- 
termined by first fitting a Gaussian profile to the Ha feature 
to obtain the redshift and FWHM, and using these values to 
constrain the fit to the [N 11] emission line. This method is 
based on the assumption that the Ha and [N 11] lines have ex- 
actly the same redshift and FWHM, with the Ha line having 
a higher signal-to-noise ratio (S/N). For most of the objects in 
our sample, the [N 11] A6548 line was too faint to measure. 
[O 111] A5007 and H/3 fluxes were determined with indepen- 
dent fits. In most cases, redshifts from Ha, [O 111] A5007, 
and H/3 agree to widiin Az = 0.0004 (Av = 50-60 km s"i at 
z = 1.0-1.4). For flie object 42021412, [O III] A5007 Ues on 
top of a bright sky line, and only a lower limit is given. 

SFRs inferred from Ha luminosities using the calibration 
of Kennicutt (1998) are shown in Table (2). The results 
have been converted from the Salpeter IMF used by Kenni- 
cutt (1998) to a Chabrier (2003) IMF by dividing flie results 
by a factor of 1 .8. In our whole sample of 20 galaxies, the Ha 
fluxes range from 5.6 x 10"'^ to 2.4 x 10"'* erg s"' cm"^. The 
mean Ha flux for the sample at z 1 .0 (z '--^ 1 .4) is 1 .3 x 10"'* 
(1.2 X 10"'*), corresponding to a star formation rate of 3 (6) 
Mq yr"', uncorrected for dust extinction or aperture effects, 
which may amount to a factor of 2 difference (Erb et al. 
2006b). Note that these characteristic Ha star formation rates, 
after being corrected for aperture effects, would be signifi- 
cantly higher than those of local galaxies in the SDSS sample 
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Fig. 2.— NIRSPEC spectra of DEEP2 galaxies in our new sample at z ~ 1.0. H/3 and [O III] are observed in the NIRSPEC 1 filter, with Ha and [N II] in the 
NIRSPEC 3 filter (similar to the J band). The 1 a error spectra are shown as dotted lines, offset vertically by -5x10"'^' 
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of Tremonti et al. (2004) and the (z) = 0.4 TKRS subsample of 
Kobulnicky & Kewley (2004), even before correction for dust 
extinction. We also note that the mean specific SFR for both 
Z - 1.0 and z - 1.4 samples is log((SFR/M^) yr"') = -9.7. 

As discussed in Paper 1, absolute line flux measurements 
suffer from several sources of systematic error, which can 
amount to at least a ~ 25% uncertainty (Erb et al. 2003). This 
level of uncertainty is present even under photometric con- 
ditions, which may not have applied through the full extent 
of our observations. For the remainder of the discussion we 
therefore focus on the measured line ratios, [N 11] A6584/Hq! 
and [O III] A5007/H/3, which are not only unaffected by un- 
certainties in flux caUbration and other systematics but also 
relatively free from the effects of dust extinction, due to the 
close wavelength spacing of the lines in each ratio. Hereafter, 
we use "[N ll]/Ha'" to refer to the measured emission-line flux 
ratio between [N II] A6584 and Ha, and "[O III]/H/?" for that 
between [O III] A5007 and H/3. 

3. THE OXYGEN ABUNDANCE 



H II region metallicity is an important probe of galaxy for- 
mation and evolution, as it represents the integrated products 
of past star formation, modulated by the inflow and outflow 
of gas. Oxygen abundance is often used as a proxy for metal- 
licity since oxygen makes up about half of the metal content 
of the interstellar medium and exhibits strong emission lines 
from multiple ionization states in the rest-frame optical that 
are easy to measure. For comparison, we use the solar oxy- 
gen abundance expressed as 12 + log(0/H) = 8.66 (Allende 
Prieto et al. 2002; Asplund et al. 2004). 

The most robust way to estimate the oxygen abundance is 
the so-called direct method, based on the measurement of 
the temperature-sensitive ratio of auroral and nebular emis- 
sion lines. However, in distant galaxies the auroral lines 
are almost always undetectable (but see Hoyos et al. 2005; 
Kakazu et al. 2007) since they become extremely weak at 
metallicities above ~0.5 solar. Even at lower metalhcities, 
the auroral lines are typically beyond the reach of the low 
S/N typical of the spectra of distant galaxies. For distant star- 
forming galaxies, therefore, measuring strong emission-line 
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ratios is the only viable way of obtaining the H II region gas- 
phase oxygen abundance (Kobulnicky et al. 1999; Pettini et al. 
2001). 

Given our NIRSPEC data set, and our desire to avoid the 
systematic uncertainties entailed in adding [O II] line fluxes 
obtained with the DEIMOS spectrograph (without real-time 
flux-calibration) to [O III] fluxes obtained with NIRSPEC, 
we focus on two strong-line ratios as indicators for the oxy- 
gen abundance: N2= log([N II]/Ha) and 03N2= log{([0 
III]/H/3)/([N II]/Ha)}. These indicators have been calibrated 
by Pettini & Pagel (2004) using local H II regions, most of 
which have direct abundance determinations. The sensi- 
tivity of the indicators to oxygen abundance, as well as their 
limitations, have been discussed in Pettini & Pagel (2004) and 
Paper I. Absolute estimates of metal abundances are quite un- 
certain, as abundances determined with different indicators 
or with different cahbrations of the same indicator may have 
substantial biases or discrepancies (e.g. Kennicutt et al. 2003; 
Kobulnicky & Kewley 2004). We therefore emphasize rela- 
tive abundances determined with the same method, using the 



same calibration. 

The N2 indicator, pointed out by several works (Storchi- 
Bergmann et al. 1994; Raimann et al. 2000; Denicolo et al. 
2002), is related to the oxygen abundance via 



12+log(0/H)= 8.90 -1-0.57 xN2, 



(1) 



which is valid for 7.50 < 12-l-log(0/H) < 8.75, with a 1 a 
scatter of ±0.18 dex (Pettini & Pagel 2004). It has been used 
by Erb et al. (2006a) to estimate oxygen abundances for UV- 
selected z ^ 2 galaxies. Note that the N2 indicator is not 
sensitive to increasing oxygen abundance above roughly so- 
lar metallicity, as shown with photoionization models (Kew- 
ley & Dopita 2002). Thus, for a subset of 12 galaxies in our 
sample with measurements of the full set of H/3, [O III], Ha, 
and [N II], we also use the 03N2 indicator introduced by Al- 
loin et al. (1979), which is expected to be particularly useful 
at solar and super-solar metalhcities where [N II] saturates 
but the strength of [O III] continues to decrease with increas- 
ing metallicity. Using the same calibration sample, Pettini & 
Pagel (2004) show that 03N2 is related to the oxygen abun- 
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dance via 

12+log(O/H) = 8.73-0.32xO3N2, (2) 
which is valid for 8.12 < 12 + log(0/H) < 9.05, with a 1 a 
scatter of ±0.14 dex. Table 2 lists oxygen abundances de- 
rived using these two indicators. The errors on the oxygen 
abundances are dominated by the systematic uncertainties in 
the calibrations of the indicators. 

3.1. Composite Spectra 

Relative, average abundances determined from composite 
spectra can be more accurately determined than those from 
individual spectra. As discussed by Erb et al. (2006a), mak- 
ing a composite spectrum not only reduces the uncertainties 
associated with the strong-hne calibration by a factor N^/'^, 
where is the number of objects included in the composite 
spectrum, but also enhances the spectrum S/N since sky lines 
generally lie at different wavelengths for spectra at different 
redshifts. In addition, one of our goals is to determine the av- 
erage properties of subgroups of galaxies in our sample. For 
the subset of 18 galaxies in our sample with Ha and [N II] 
measurements, as well as stellar mass estimates, we divide 
the sample into four bins by stellar mass, with two bins at 
z ~ 1.0 and two bins at z ~ 1.4. For the subset of 12 galaxies 
with measurements of not only Ha and |N II |, but also H/3 
and [O III], we also divide the sample into four bins by stellar 
mass with two bins each at z 1.0 and at z ~ 1 .4. 

To make the composite spectra, we first shift the individ- 
ual one-dimensional flux-calibrated spectra into the rest frame 
and then combine them by generating the median spectrum, 
which preserves the relative fluxes of the emission features 
(Vanden Berk et al. 2001). We use N2 and N2+03 com- 
posite spectra to refer to the composites with Ha and [N 
II], and those with all four lines, respectively. The N2 and 
N2+03 composite spectra, labelled with mean stellar mass in 
each bin, are shown in Figures 4 and 5. The corresponding 
emission-line flux ratios along with uncertainties measured 
from the composite spectra, as well as the inferred oxygen 
abundances, are hsted in Tables 3 and 4. The listed errors in 
12 + log(0/H) include the uncertainties from the propagation 
of emission-line flux measurements, as well as the systematic 
scatter from the strong-line calibration. As shown in paper 
I, the systematic discrepancies between the N2- and 03N2- 
based abundances are mainly due to the fact that, on aver- 
age, DEEP2 galaxies are offset from the excitation sequence 
formed by local H II regions and star-forming galaxies. We 
discuss this issue in detail in §5. 

4. THE MASS-METALLICITY RELATION 

The redshift evolution of the luminosity-metallicity and 
mass-metalUcity relations provides important constraints on 
models of galaxy evolution. A correlation between gas-phase 
metallicity and stellar mass can be explained by either the ten- 
dency of lower mass galaxies to have larger gas fractions and 
lower star formation efficiencies (McGaugh & de Blok 1997; 
Bell & de Jong 2000; Kobuhiicky et al. 2003), or the preferen- 
tial loss of metals from galaxies with shallow potential wells 
by galactic-scale winds (Larson 1974). In the local universe, 
strong correlations between rest-frame optical luminosity and 
the degree of chemical enrichment have been observed in both 
star-forming and early-type galaxies (Garnett & Shields 1987; 
Brodie & Huchra 1991; Tremonti et al. 2004). The correla- 
tion has also been observed in intermediate- and high- red- 
shift samples (Kobulnicky et al. 2003; Lilly et al. 2003; Kob- 
ulnicky & Kewley 2004; Erb et al. 2006a), altiiough caution 



must be taken when comparing samples with metallicities de- 
termined from different methods. Physically, the correlation 
between stellar mass and metallicity is more fundamental than 
that between luminosity and metaUicity (Paper I; Tremonti 
et al. 2004; Erb et al. 2006a). We therefore focus on the mass- 
metallicity relation in the following discussion. 

The left panel of Figure 6 shows the average metallicity 
of the galaxies in each mass bin determined from the N2 
composite spectra plotted against their average stellar mass 
at z ^ 1.0 {open diamonds) and at z ^ 1-4 {open squares). 
Although our sample is still small, we do see evidence for 
mass-metallicity relations at both z ^ 1.0 and z ^ 1.4. These 
trends are also present when we examine the metallicities and 
stellar masses for individual objects, which are plotted in the 
figure as weU. For comparison, the local SDSS galaxies dis- 
cussed by Tremonti et al. (2004) are denoted by contours and 
dots^ and the z 2 Erb et al. (2006a) sample as open cir- 
cles. Metallicities for SDSS galaxies were calculated using 
the same strong-line indicator that was applied to the DEEP2 
galaxies and not the Bayesian O/H estimate from Tremonti 
et al. (2004). 

At fixed stellar mass, the metallicities of our z ^ 1.0-1.5 
sample as a whole are lower than those of local galaxies yet 
higher than those of the z ~ 2 sample. However, there is evi- 
dence for a reverse trend between the subsets of our sample at 
z ~ 1 .0 and at z ~ 1 -4. In Paper I, this difference was attributed 
to the fact that the z 1 .0 sample was on average fainter and 
less massive than the z ^ 1 .4 sample. With a larger sample, 
however, we find that the higher mass z 1 bin does have 
lower metallicity than the lower mass z ~ 1-4 bin. Differ- 
ences in outflow or inflow rate of unenriched gas at z ~ 1-0 
and at z ^ 1.4 could give rise to this trend. However, the in- 
terval in cosmic time between z ~ 1 4 and 1 .0 is small enough 
that typical gas inflow rates at fixed stellar mass, and the cor- 
responding star formation and outflow rates, will not signifi- 
cantly evolve. Therefore, this explanation is not a Ukely cause 
of the reverse trend. A different average degree of dust red- 
dening at z ~ 1-0 and at z 1 .4 is also not a likely cause, 
since the N2 indicator is based on emission lines with very 
close spacing in wavelength. On the other hand, if the z ~ 1 .0 
galaxies have systematically different physical conditions or 
less significant contributions from AGN activity relative to the 
z ^ 1 .4 objects, metallicities estimated with the same calibra- 
tion would be systematically biased between the two samples 
in such a way to produce the observed trend. As discussed in 
§5, we propose that the most likely cause for the reverse trend 
is this difference in H II region physical conditions. Since the 
systematic uncertainty from strong-line calibration is large, 
and our sample is still too small to draw any solid conclu- 
sion, it will become feasible to clarify this issue only when a 
statistically large enough sample is assembled, and both the 
high- and low- mass ends are spanned at z ^ 1.0 as well as at 
z- 1.4. 

We also plot the mass-metallicity relation from the N2H-03 
composite spectra in Figure 6, where the left panel shows 
metallicities determined from the N2 indicator, while the right 
panel shows those determined from the 03N2 indicator, at 
bothz~ I. Q {filled diamonds) andz-^ 1.4 {filled squares). The 
03N2-based abundances are systematically lower than those 

^ Note that here and throughout, we use a combination of contours and 
dots to indicate SDSS objects. On each such appHcable plot, SDSS data 
points were mapped onto 10 evenly spaced levels according to surface den- 
sity, where objects on the lowest level are denoted by dots while other levels 
are presented by contours. 
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Fig. 4. — Composite NIRSPEC spectra of the 18 z ~ 1 .0- 1 .5 galaxies with both stellar mass estimation as well as [N II] and Ha measurement in our sample, 
divided separately at 7, ~ 1.0 and 7. ^ 1.4. The 1 cr error spectra are shown as dotted lines, offset vertically by -3 X Ur"* ergs s"' cm"^ A"' for clarity. The 
spectra are labeled with the mean stellar mass from each bin, and the Ha, [N II] and [S II] lines are marked by dashed lines. Note that the spectra near the 
density-sensitive [S 11] lines are very noisy, due to the large dispersion of the flux per count from near-IR standard star calibration, caused by low efficiency near 
the filter edge. 



based on N2. As suggested in Paper I and discussed in detail 
in §5, these systematic discrepancies between N2- and 03N2- 
based abundances are due to the fact that DEEP2 galaxies de- 
part from the local H 11 region excitation sequence. In ad- 
dition, the reverse trend between z ~ 1 .0 and 1 .4 in metal- 
licity estimated from 03N2 is much less significant than the 
one in metallicity estimated from N2. We return to this issue 
as well in §5. Despite these discrepancies, the overall corre- 
lation between average stellar mass and metallicity observed 
among the N2 composite spectra is still present for the N2+03 
spectra. This is evidence that the correlation is insensitive to 
the spectrum of any particular object, as it is robust to anal- 
yses using different binning schemes. At the lower mass end 
(M,, ~ 8 X 10'' Mq), the average metallicity of z ~ 1.0-1.5 
galaxies based on N2 is at least 0.22 dex lower than the lo- 
cal typical value. Since the N2 indicator saturates near solar 
abundance, as discussed in Erb et al. (2006a), it is difficult to 
determine the true metallicity offset between two samples at 
different redshifts using this indicator. We can also determine 
the offset from the 03N2-based abundances, particularly near 
the solar abundance. Based on 03N2 abundances, the metal- 
licity offset between our DEEP2 objects and the local SDSS 
sample is at least 0.21 dex at the high-mass end {M^, ~ 5 x lO"* 
Mq). However, as we discuss in §5, these estimates based on 
strong-line indicators may still be subject to systematic uncer- 
tainties from using the calibration of H n regions with signif- 
icantly different physical properties. 

5. THE OFFSET IN DIAGNOSTIC LINE RATIOS OF HIGH-REDSHUT 
GALAXIES 

5.1. Emission-Line Diagnostics 

There is evidence that physical conditions in the H II re- 
gions of high-redshift galaxies hosting intense star forma- 
tion are different from those of local H II regions (Paper I). 



TABLE 3 

Oxygen Abundances from N2 composite spectra. 



Bin 


N" 




log(M, /M^f 


N2'' 


12 + logCO/Hf 


1 

2 

3 

4 


3 
4 
5 
6 


1.0375 
1.0207 
1.3930 
1.3903 


9.88:0;|« 
10.37:!];«i 
lO.lO^o* 
10.83^ 


-0.96±0.03 
-0.74±0.02 
-0.63±0.02 
-0.49±0.02 


8.35±0.11 
8.48±0.09 
8.54±0.09 
8.62±0.08 



^Number of objects contained in each bin. 
''Mean redshift for each bin. 

'^Mean stellar mass and uncertainty from error propagation. 
''N2 = log([N n]A6584/Ha). 

'^Oxygen abundance deduced from the N2 relationship presented in Pettini & Pagel 
(2004). 

The most common method for probing H II region physical 
conditions, and discriminating between star-forming galax- 
ies and AGNs, is based on the positions of objects in the 
Baldwin, PhilUps, & Terlevich (1981) empirical diagnostic 
diagrams (hereafter BPT diagrams). These plots feature the 
optical line ratios fN Ill/Ha, [O I]/Ha, [S III/Hq, and [O 
IlI]/H/3, and have been updated by many authors, including 
Osterbrock & Pogge (1985), Veilleux & Osterbrock (1987), 
Kewley et al. (2001a, hereafter KeOl), and Kauffmaim et al. 
(2003a, hereafter Ka03). A considerable fraction of high- 
redshift star-forming galaxies at both z ^ 1 (paper I) and z 
> 2 (Erb et al. 2006a) do not follow the local excitation se- 
quence formed by nearby H 11 regions and star-forming galax- 
ies on the emission-line diagnostic diagram of [N II]/Ha ver- 
sus [O in]/H/3; on average, they lie offset upward and to the 
right. These differences must be taken into account when ap- 
plying empirically calibrated abundance indicators to galaxy 
samples at different redshifts. Several possible causes for the 
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Fig. 5. — Composite NIRSPEC spectra of the 12 z ~ 1 — 1.5 galaxies with both stellar mass estimation as well as all-four-line measurement in our sample, 
divided separately at z ~ 1.0 and z ~ 1.4. The 1 <t error spectra are shown as dotted hnes, offset vertically by -3 x 10"'* ergs s"' cm"^ A~' for clarity. The 
spectra are labeled with the mean stellar mass from each bin, and the H/3, [O HI], Ha, [N 11] and [S n] lines are marked by dashed Unes. 



offset have been discussed in Paper I, including differences in 
the ionizing spectrum, ionization parameter, electron density, 
and the effects of shocks and AGNs. 

In Figure 7, fO llll/H/3 and [N 11]/Hq! ratios are plotted in 
the left panel for the 13 galaxies in our sample with the full set 
of emission lines, and the average [O lll]/H/3 and [N 11]/Hq! 
ratios from the N2H-03 composite spectra are also shown in 
the right panel for bins at z ~ 1 .0 and 1 .4, respectively. A 
subset of emission-line objects from SDSS are also shown as 
grey contours and dots for comparison. The dotted curve is 
from KeOl, representing a theoretical upper limit on the loca- 
tion of star-forming galaxies in the diagnostic diagram. The 



dashed curve is from Ka03 and serves as an empirical discrim- 
inator between star-forming galaxies and AGNs. On average, 
the Ha flux of galaxies below this curve should have < 1% 
contribution from AGNs (Brinchmann et al. 2004). The effect 
observed in Paper 1 is still present in our larger sample, in the 
sense that the z ^ 1 .0- 1 .4 sample is, on average, significantly 
offset from the excitation sequence formed by star-forming 
galaxies from SDSS. Furthermore, the average offset for the 
z 1.4 objects is larger than for those at z 1.0. A similar, 
if not stronger, effect is observed in star-forming galaxies at 
z ^ 2 (Erb et al. 2006a). As discussed in Paper 1, unaccounted- 
for stellar Balmer absorption is not the explanation for the off- 
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log (M./Ms) log (M./Mq) 

Fig. 6. — Mass-metallicity relation observed at z ~ 1.0 and z ~ 1-4. In the left panel, metallicities are infeiTed from the N2 indicator, while the right panel 
shows metallicities estimated from the 03N2 indicator In both plots, open (filled) diamonds and squares are used for N2 (N2+03) composite spectra at z ~ 1.0 
and 1.4, respectively (see section 3.1). EiTors show the uncertainties propagated from the strong-line ratio measurement, while the systematic uncertainties 
from the strong-Hne method calibration (Pettini & Pagel 2004) are shown in the lower right comer of each plot as dotted (solid) error bars for N2 (N2-I-03) 
composite spectra. The systematic uncertainties from the calibration are reduced by a factor W'/^, where N is the number of individual spectra included in each 
composite. For stellar masses, the horizontal bars give the mass range in each bin. Data points for individual objects are shown as filled triangles and circles for 
z ~ 1.0 and z ~ 1-4, respectively. Associated error bars are listed in Table 2 in both this paper and in Paper I. For comparison, metallicities as a function of stellar 
mass are also shown, for both local SDSS galaxies (grey contours and dots) and the Erb et al. (2006a) z ~ 2 sample {grey open circles, left panel). Note that here 
and throughout, we use contours and dots to show SDSS objects. On each such applicable plot, SDSS data points were mapped onto 10 evenly spaced levels 
according to number density, where objects on the lowest level are denoted by dots while other levels are presented by contours. Solar metallicity is marked with 
a horizontal dashed line. It can be seen that the N2 indicator saturates near the solar abundance. 



TABLE 4 

Oxygen Abundances from N2-I-03 composite spectra. 



Bin N (zHa) log(M,/MQ) N2 03N2" [12 + log(0/H)]jv2'' [12 + log(0/H)]o3iV2'-' 



1 3 1.0372 9.98:™09 -0.89±0.03 1.37±0.03 8.39±0.10 8.29±0.08 

2 3 1.0216 10.40;g||^ -0.79±0.02 1.17±0.06 8.45±0.10 8.35±0.08 

3 3 1.3923 ^"•''6-006 -0.61±0.04 1.09±0.08 8.55±0.11 8.38±0.08 

4 3 1.3858 10.63!™^ -0.51±0.02 > 0.62 8.61±0.10 < 8.53 



"03N2 = log{([0 m]A5()07/H/3)/([N II] A6584/Hq)}. 

^Oxygen abundance deduced from the N2 relationship presented in Pettini & Pagel (2004). 
'' Oxygen abundance deduced from the 03N2 relationship presented in Pettini & Pagel (2004). 



sets in emission-line ratios, since the corrections would shift 
the DEEP2 galaxies by no more than 0. 1 dex downward and 
by an insignificant amount in [N II]/Ha. 

Isolating the causes of the offset of the high-redshift sam- 
ples from the local excitation sequence on the diagnostic di- 
agram will provide important insight into the physical condi- 
tions in distant star-forming galaxies. These conditions also 
comprise an essential systematic factor determining the emer- 
gent strong emission-line ratios and associated calibration of 
chemical abundances. Unfortunately, we do not currently 
have a large high-redshift sample with available additional 
spectral features, stellar population parameters, and morpho- 
logical information, which is needed for a direct study of how 
the diagnostic line ratios vary as functions of galaxy proper- 
ties. However, as seen in Figure 7, while only a tiny fraction 
of the SDSS galaxies reside in the region of BPT parame- 
ter space inhabited by our most extremely offset high-redshift 
galaxies, in between the star-forming excitation sequence and 
the AGN branch, the sheer size of the SDSS sample still re- 
sults in a set of ^ 100 such local objects. These objects 
can serve as possible local counterparts for our DEEP2 ob- 
jects, with the added benefit of high S/N photometric, spectro- 



scopic, and morphological information from SDSS. We will 
make use of this detailed information to understand the cause 
of the local galaxies' offset in the BPT diagram, and, by ex- 
tension, the likely cause of the offset among the high-redshift 
galaxies. 

In the following sections, we compare in detail these 
anomalous SDSS objects against more typical SDSS star- 
forming galaxies. Accordingly, we analyze possible causes 
for their offset in terms of different physical conditions in the 
ionized regions, which include H 11 region electron density, 
hardness of the ionizing spectrum, ionization parameter, the 
effects of shock excitation, and contributions from an AGN. 
We further try to unravel possible connections between physi- 
cal conditions of these anomalous SDSS objects and their host 
galaxy properties and use them to interpret our observations 
of high-redshift star-forming galaxies. 

5.2. Local Counterparts: SDSS Main and Offset Samples 

Within this work we use the SDSS Data Release 4 (DR4; 
Adelman-McCarthy et al. 2006) spectroscopic galaxy sam- 
ple, which contains m-, g-, r-, /-, and z-band photometry and 
spectroscopy of 567,486 objects. As described in Tremonti 
et al. (2004), emission-line fluxes of these galaxies were mea- 



STAR-FORMING GALAXIES AT Z~ 1 .0- 1 .5 11 



to. 



o 
o 



1 .5 



1.0 



0.5 



0.0 



-0.5 



1 1 1 1 1 1 1 
" o pilot, z ~ 1 .0 


1 1 1 1 1 1 1 1 1 1 1 


- □ pilot, z ~ 1 .4 




" ♦ new, z ~ 1 .0 


- 


-■new, z ~ 1 .4 
















. . . . 1 . . 


- I 

. 1 . . 1 . 1 . . . 



to. 



o 
o 



1 .5 



1.0 



0.5 



0.0 



-0.5 



-1.0 



▲ K 1 n c n t 
▼ u 1 1 1 U I 


1 1 1 1 1 
z ~ 1 


1 1 1 1 1 1 1 1 


- ■ bins at 


z ~ 1 .4 


- 


- 




- 




1 


\ \\ \\ 

' hi Im::.^. -■- 










.... 1 


\ \ 
,1,1,,.,, 



-1.5 -1.0 -0.5 0.0 
log ([Nll]6584/Ha) 



0.5 



■1.5 -1.0 -0.5 0.0 
log ([Nll]6584/Ha) 



0.5 



Fig. 7. — H n region log[N II]/Ha vs. log[0 III]/H/3 diagram. On tlie left, tlie z ~ 1.0- 1.5 DEEP2 galaxies are shown as indicated on the plot. The right panel 
shows the average diagnostic line ratios from the composite spectra of DEEP2 objects. Note the bin with the highest [N II]/Ha value only has a lower limit on [O 
III]/H/3 because it includes the object 42021412. SDSS emission-line galaxies that satisfy the criteria described in §5.2 ai'e shown as grey contours and dots. The 
dashed line is an empirical demarcation from Ka03 based on the SDSS galaxies, whereas the dotted line is the theoretical limit for star-forming galaxies from 
KeOl. Nearby star-forming galaxies and H II regions form a well-defined excitation sequence of photoionization by massive stars, below and to the left of these 
curves. The z ~ 1.0—1.5 DEEP2 objects on average are offset from this excitation sequence, with objects at z ~1.4 more offset than those at z ~1.0. 

sured from the stellar-continuum subtracted spectra with the 
latest high spectral resolution population synthesis models by 
Bruzual & Chariot (2003). 

Before being further divided into groups with different 
emission-line diagnostic ratios, our sample was selected from 
the 567,486 galaxy DR4 sample according to the following 
criteria: 

(1) Redshifts in the range 0.005 <z< 0.25. (2) Signal- 
to-noise ratio (S/N) > 3 in the strong emission-lines [O II] 
AA3726, 3729, H/3, [O III] AA4959, 5007, Ha, [N II] A6584, 
and [S II] AA6717, 6731; and S/N > 1 in the weak emission- 
line [O I] A6300. In practice, ^ 99% of the objects also have 
[O I] A6300 S/N>3. (3) The fiber aperture covers at least 20% 
of the total ,g-band photons. (4) Stellar population parameter 
estimates are available in the catalog derived using methods 
described by Kauffmann et al. (2003b). 

The first criterion on redshift range is the same as that 
adopted by Tremonti et al. (2004), which enables a fair 
comparison with the mass-metallicity relation from their 
work. The second criterion selects galaxies with well- 
measured emission-line fluxes, required by our study based 
on emission-line diagnostic ratios. The S/N values were ob- 
tained using errors of the emission-line fluxes, which were 
first taken from the emission-line flux catalog on the MPA 
DR4 Web site^ and then scaled by the recommended factors 
(see the Web site for more details about the emission-line 
flux error). The third criterion is required to avoid signifi- 
cant aperture effects on the flux ratios (Kewley et al. 2005; 
Tremonti et al. 2004). Lower aperture fraction can cause sig- 
nificant discrepancies between aperture and global parameter 
estimates. The fourth criterion enables a comparison of stel- 
lar population properties with diagnostic emission-line ratios. 
Roughly 30% of the initial set of DR4 galaxies were ruled 
out because they do not have stellar population parameter es- 
timates; an additional ^20% was cut out according to the red- 



' See http://www.mpa-gai'ching. mpg.de/SDSS/ 



shift and aperture-coverage constraints. Rules on the S/N in 
both strong and weak emission lines removed another ^45% 
of the objects. In all, these selection criteria leave us with 
^31,000 emission-line objects, shown as grey contours and 
dots in Figure 7, which is about 5% of the 567,486 galaxy 
DR4 sample. 

With the exception of one of the 13 DEEP2 objects, which 
lies on the part of the [N II]/Hq; versus [O III]/H/3 diagram 
populated by many H II/AGN composites, the high-redshift 
galaxies on average fall between the excitation sequence of 
typical SDSS star-forming galaxies and the AGN branch. Al- 
though many of the DEEP2 objects have a less extreme offset 
relative to typical SDSS star-forming galaxies than those of 
the sample of Erb et al. (2006a), and the most extreme z ^ 1 .4 
objects in Paper I, there are still several that are significantly 
offset from the main sequence. We proceed by studying prop- 
erties of the SDSS galaxies with similar [N II]/Ha and [O 
III]/H/3 values to those of the significantly offset DEEP2 ob- 
jects, and comparing these unusual objects with typical SDSS 
star-forming galaxies along the excitation sequence. In this 
sense we further divide the SDSS sample into the "Main", 
and "Offset" subgroups, as shown in Figure 8. 

Objects in the Main sample are selected to lie below the 
Ka03 empirical curve, whereas Offset objects are located in 
between the Ka03 and the KeOl curves. The Offset sample 
is also constrained to have log([N II]/Hq;) < -0.44 and log([0 
III]/H/3) > 0.26, according to the two emission-line ratios of 
the DEEP2 object with the second largest [N II]/Ha value. 
We do not use the galaxy with the largest [N II]/Ha value in 
our DEEP2 sample for the selection criteria of the two line 
ratios, because the number density in the composite region 
increases very rapidly as the AGN branch is approached, yet 
most of our DEEP2 galaxies clearly do not fall in that regime. 
Also, it is worth pointing out that, while they are significantly 
displaced from the local excitation sequence, more than half 
of the DEEP2 objects actually lie below the Ka03 curve. In- 
deed, the SDSS Offset sample is constructed according to our 
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Fig. 8. — H II region diagnostic diagram: log[N II]/Ha vs. log[0 III]/H/3. SDSS Main, Offset- AGN, Offset-ambiguous, and Offset-SF samples are shown as 
grey contours with dots, blue triangles, green crosses, and red circles, respectively (see §5.2 for more information on sample selection rules). The z ~ 1.0- 1.5 
DEEP2 objects are also shown for a comparison using the same legends as in the left panel of Figure 7, and the dotted and dashed lines have the same meanings 
as those in Figure 7. 



DEEP2 objects that are offset the most, in order to create a 
strong contrast with the Main sample. It is also representative 
of the z ^ 2 galaxies observed by Erb et al. (2006a). As we de- 
scribe in §5.6, the conclusions drawn from this extreme sam- 
ple are corroborated by the work of Brinchmann et al. (2008), 
where objects with less extreme offsets are considered and 
therefore should be valid for typical objects in our DEEP2 
sample. Finally, the Offset sample only covers a certain range 
of stellar masses, so we further select the Main control sample 
according to the same stellar-mass range in order to ensure a 
fair comparison. This leaves us with ^^2 1,000 objects for the 
Main sample and 101 objects for the Offset sample. 

Figure 9 displays how objects in the Offset sample are dis- 
tributed in the additional BPT diagrams featuring [O I]/Hq! 
and [S II]/Ha. These plots indicate that Offset objects se- 
lected solely on the basis of their [O III]/H/3 and [N II]/Ha 
ratios span a diverse range of properties in other physical pa- 



rameter spaces. High [O I]/Ha and [S II]/Hq; both occur 
when there is a hard ionizing radiation field, including sig- 
nificant contribution from X-ray photons. In this case, there 
is an extended, partially ionized zone, where H I and H II 
coexist, and [O I] and [S II] are dominant forms of O and 
S. The extended zone of partially ionized H does not exist 
in H II regions photoionized by OB stars (Evans & Dopita 
1985; Veilleux & Osterbrock 1987). High [O I]/Ha and [S 
II]/Hq! are also produced in gas that has been heated by fast, 
radiative shocks, which also produce partially ionized shock- 
precursor regions (Dopita & Sutherland 1995, 1996). Mate- 
rial in supernova remnants provides an example of shocked 
gas. While it is sensitive to shocks, [S II] is more suscep- 
tible to collisional de-excitation than [O I], given its critical 
density (2x10^ cm"^), as opposed to that of [O I] (2x10'' 
cm"-'), and therefore might be suppressed in regions of high 
electron densities (Dopita 1997; Kewley et al. 2001b). An- 
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Fig. 9. — H II region diagnostic diagrams. Descriptions for SDSS objects are tlie same as those in Figure 8. On tlie [O I]/Hq and [S II]/Ha diagrams, tlie dotted 
(blue solid) lines are empirical curves separating star-forming galaxies (Seyferts) and AGN (LINERs) from Kewley et al. (2006). The log[N II]/[0 II] vs. log[0 
III]/H/3 diagram is useful in separating ionization parameter and abundance but becomes insensitive to ionization parameter above solar metalhcity (Dopita et al. 
2000). The [N II]/[0 II] flux ratio has been corrected for dust extinction, given the wide spacing in wavelength between [N II] and [O II]. 



Other point worth mentioning is that [O IJ/Ha should reveal 
larger differences in regions ionized by hard spectrum as op- 
posed to that of stars, than [S II]/Hq; does. This is because 
the ionization potential of [O I] matches that of H I, so [O 
I] enhancement should reflect increased presence of partially 
ionized zone. However, [S II] exists in both completely and 
partially ionized zones. So the contrast between AGN- and 
stellar-ionized regions for [S II]/Hq; is not as great as that for 
[O I]/Ha. 

In summary, the various locations of the Offset objects on 
the [O I]/Ha and [S II]/Hq; diagrams indicate different levels 
of contribution from AGN and shock excitations to the emerg- 
ing spectra. We therefore further divide the Offset sample on 
the [O l]/Ha and [S II]/Hq; diagrams, according to the the- 
oretical scheme for classifying starburst galaxies and AGNs 
(Kewley et al. 2001a, 2006). It is worth noting here that all ob- 
jects in the Offset sample have S/N>3 in both [S II] and [O I] 
emission lines, so this division should not be compromised by 
spurious measurements. While there may still be a low-level 
(< 10%) AGN contribution to the spectra of objects classi- 
fied as starbursts by this scheme, it serves as a rough guide to 
the range of properties in the Offset sample. This classifica- 
tion leaves us with (1) "Offset-AGN" (43), where objects lie 
above both the KeOl curves of the two diagrams; (2) "Offset- 
ambiguous" (33), where objects lie above one KeOl curve and 
under the other one; and (3) "Offset-SF" (25), where objects 
lie below both the KeOl curves. The division of Offset-AGN 
and Offset-SF is only based on hardness of the ionizing spec- 
trum and the contribution of shock excitation. We show in 
Section 5.4 that these two offset samples have different host 
galaxy properties. 

5.3. H II Region Emission-Line Diagnostic Ratio, Physical 
Conditions, and Galaxy Properties 

In order to determine the origin of the observed difference 
in emission-line ratios between the IVIain and Offset samples, 
we consider a large set of galaxy properties. Measurements 
of emission lines as well as photometric, physical, and envi- 
ronmental properties of stellar population are available for our 
SDSS samples. In this section, we describe the galaxy prop- 
erties of interest, while in the following section, a comparison 
of the Main and Offset sample distributions in these properties 
is presented. 

First, using emission-line ratios, we examine the main fac- 
tors controlling the emission-line spectrum in an H II region, 
i.e., the gas phase metallicity, the shape or hardness of the ion- 



izing radiation spectrum, and the geometrical distribution of 
gas with respect to the ionizing sources, which can be repre- 
sented by the mean ionization parameter and electron density 
(Dopita et al. 2000). Several H II region emission-line ratios 
can be used as diagnostics of these factors. 

Apart from the N2 and 03N2 indicators, we also use the 
ratio N202 = log([N II] A6584/[0 II] AA3726, 3729), sug- 
gested by van Zee et al. (1998), as an abundance diagnos- 
tic for SDSS objects, because N202 is virtually independent 
of ionization parameter and also strongly sensitive to metal- 
licity. This indicator is monotonic between 0.1 and over 3.0 
times solar metallicity (Dopita et al. 2000; Kewley & Dopita 
2002). Concerns about reddening correction and reliable cali- 
bration over such a large wavelength baseline have hampered 
the use of this ratio. As shown in Kewley & Dopita (2002), 
however, the use of classical reddening curves and standard 
calibration are quite sufficient to allow this [N II]/[0 II] di- 
agnostic to be used as a reliable abundance indicator The 
Bresolin (2007) calibration of N202 is used to infer oxygen 
abundance. We also use the R2i = log{([0 II] A3727 + [O 
III] AA4959, 5007)/H/3} pai-ameter, introduced by Pagel et al. 
(1979), as an abundance indicator, to make comparison with 
several works (Lilly et al. 2003; Savaglio et al. 2005), al- 
though this indicator has some well-documented drawbacks 
(e.g. Kobulnicky et al. 1999; Kewley & Dopita 2002). 

C»32 = log([0 III] A5007/[O II] AA3726, 3729) is used as 
an indicator of the ionization parameter (McGaugh 1 99 1 ; Do- 
pita et al. 2000). However, is sensitive to both ioniza- 
tion parameter and abundance; both low ionization parameter 
and high abundances produce low values of 032 (Dopita et al. 
2000; Kewley & Dopita 2002). Thus, when 032 is used as 
a diagnostic of the relative ionization parameters of different 
samples, the metallicities of these samples must also be taken 
into account for a fair comparison. A diagnostic plot, for ex- 
ample, N202 versus 032, can be used to separate ionization 
parameter and abundance. The ratio [S III]/[S II] provides 
another diagnostic of the ionization parameter, which is inde- 
pendent of metallicity except at very high values of the ratio 
(Kewley & Dopita 2002). Unfortunately, the SDSS spectra 
do not cover [S III] AA9069, 9532 from SDSS, so it is not 
included here. 

FinaUy, [S II] A6717/[S II] A6731 is adopted as an electron- 
density indicator We do not use the density-sensitive ratio, [O 
II] A3726/[0 II] A3729, since this doublet is typically blended 
in SDSS spectra. 

We note that emission-line ratios, such as 032, N202, and 
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Offset-SF samples (red sotid tine; see §5.2 for more information). For clarity, the distribution functions are normalized to three different levels, with the maxima 
being 1.0, 0.7, and 0.5 for the SDSS Main, Offset-AGN, and Offset-SF samples, respectively. 



R23 have large separations in wavelength and thus need to be 
corrected for dust reddening. The Balmer decrement method 
with the reddening curve of Calzetti et al. (2000) is used to 
correct these quantities for dust extinction, assuming an in- 
trinsic ratio of Ha/Hf] = 2.78, appropriate for = 10,000 
K. 

Next, we turn to the question of the stellar populations and 
structural and environmental properties of the Main and the 
Offset sample host galaxies. To examine stellar ages, we 
adopt the narrow definition of the 4000 break denoted as 
D„(4000) (Balogh et al. 1999), which is small for young stel- 
lar populations and large for old, metal-rich galaxies (Kauff- 
mann et al. 2003b). We use the HSa index as an indicator 
of recent starburst activities (e.g. Kauffmann et al. 2003b; 
Worthey & Ottaviani 1997). Strong H6 absorption lines 



arise in galaxies that experienced a burst of star formation 
that ended ^0.1-1 Gyr ago. For galaxy morphology, we 
use the standard concentration parameter defined as the ra- 
tio C = Riioir) /Rsoif), where R^oir) and Rsoir) are the radii 
enclosing 90% and 50% of the Petrosian r-band luminos- 
ity of the galaxy. For galaxy size, we study Rsoiz), the ra- 
dius enclosing 50% of the Petrosian z-band luminosity of a 
galaxy. We also examine the surface mass density /i^,, de- 
fined as M^/llTrRloiz)] (Kauffmann et al. 2003c), the SFR 
surface density 'S^ = SFR/[27r/?5Q(r)], where we use a radius 
defined in the r band rather than the z band, as it is more 
appropriate for Ha luminosities. Additional parameters in- 
clude the fiber-corrected specific star formation rate SFR/M*, 
Ha equivalent width (EW), and the color-magnitude diagram 
(g_/)0.i versus Mr, where (g-/)"' denotes the (g-i) color 
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^-corrected to z = 0.1, and Mr stands for the fc-corrected r- 
band absolute magnitude (Kauffmann et al. 2003b). In or- 
der to search for any environmental dependence of the Offset 
sample properties, we count the number of spectroscopicaUy 
observed galaxies that are located within 2 Mpc in projected 
radius and ±500 km s"' in velocity difference from each ob- 
ject in our SDSS samples. Here we follow the procedures 
described in Kauffmann et al. (2004), constructing a volume- 
limited tracer sample using the SDSS DR4 data. 

5.4. Comparison with Typical SDSS Star-forming Galaxies 

In this section, we compare the properties of Offset and 
Main sample objects, in terms of their H II region physical 
conditions, stellar populations, structural parameters, and en- 
vironments. In particular, we find striking differences in H II 
region ionization parameter, electron density, galaxy size, and 
star formation rate surface density. Figure 10 shows rela- 
tive distributions of diagnostic-line ratios and galaxy proper- 
ties of the SDSS Main, Offset- AGN, and Offset-SF samples. 
For clarity, the distribution functions are normalized to three 
different levels, with the maxima being 1 .0, 0.7 and 0.5 for 
the SDSS Main, Offset-AGN, and Offset-SF samples, respec- 
tively. 

First, we consider H II region physical conditions. Com- 
pared to the Main sample, both Offset samples on average 
have larger ionization parameters. The Offset-SF and Offset- 
AGN samples have O32 = 0.0 and -0.2, respectively, as op- 
posed to -0.7 for the Main sample. Here and throughout, we 
refer to the median values of the distributions. Since O32 also 
depends on metaUicity, we need to make a fair comparison 
for the ionization parameter at a fixed metaUicity. Figure 1 1 
shows 032 as a function of several emission-line diagnostic 
ratios and galaxy properties. In the lower four panels, we can 
see that, regardless of which metalUcity indicator is used, the 
Offset-SF sample has bigger O32 values than the Main sample 
at a fixed strong-line indicator value. Therefore, the Offset- 
SF sample has larger average ionization parameter than Main 
sample objects having similar metallicities, independent of 
the metaUicity indicator. Compared with the Main sample, aU 
indicators except R23 suggest that the Offset-AGN has higher 
ionization parameters at fixed metalUcities. 

The Offset-SF sample also has significantly larger electron 
densities, on average, than the Main sample ([S II] A6717/[S 
II] A6731 of ~1.23 compared to ~1.40, corresponding to 
electron densities of ~208 cm"^ compared to ~47 cm~^), 
whereas the Offset-AGN sample has only sUghtly higher elec- 
tron densities than those of the Main sample ([S II] A6717/[S 
II] A6731 of ~1.37, corresponding to electron densities of 
~67 cm~^). As discussed in Paper I, for fixed ionization pa- 
rameter, metalUcity, and input ionizing spectrum, the pho- 
toionization models presented in Kewley et al. (2001a) dis- 
play a dependence on electron density, in the sense that model 
grids with higher electron density have an upper envelope in 
the space of [O III]/H/3 versus |N II]/Hq! that is offset upward 
and to the right, relative to model grids with lower electron 
density. This theoretical shift in the BPT diagram due to in- 
creased electron density is quaUtatively reflected in the prop- 
erties of Offset-SF objects. 

Next, we consider photometric and spectroscopic stellar 
population properties. The Offset-SF sample is similar to 
the Main sample in terms of median color [(g - j')" ' ~ 
0.64 compared to ^0.66], stellar age [the same value of 
D„(4000) = 1.24], and fiber-corrected specific star forma- 
tion rates (the same value of log{{SFR/M^,) yr~^) = -9.7), 



whereas the Offset-AGN sample has redder colors [(g-i) of 
--0.86], older steUar ages [D„(4000) = 1.50], and lower spe- 
cific star formation rates (log((5F/?/M^) yr"') = -10.0). Also, 
the Offset-SF sample has a larger Ha EW than the Main sam- 
ple (62 compared to 40 ), whereas the Offset-AGN sample 
has a much smaller Ha EW (15). The Offset-SF sample has 
a smaller burst fraction in stellar mass than the Main sample 
(HSa of 4.3 compared to 5.5), while the discrepancy between 
the Offset-AGN and Main samples is even larger (H^a of 3.3 
compared to 5.5). 

In terms of galaxy structure, the Offset-SF sample has 
larger concentration (C = 2.58 compared to 2.36), smaller 
half-light radii [/?5()(z) of 1 .2 kpc compared to 1 .9 kpc], higher 
surface stellar mass density [log(/i*) of 8.78 Mq kpc"^ com- 
pared to 8.61 M0 kpc"^], and higher SFR surface density 
[log(S*) of -0.86 Mq yr"' kpc"^ compared to -1.12 Mq yr"' 
kpc"^] than the Main sample, whereas the Offset-AGN has 
larger concentration (C of 2.66), moderately smaller sizes 
[^5o(z) of 1.5 kpc], moderately higher surface steUar mass 
density [log(/i^) of 8.68 Mq kpc~^], and similar SFR surface 
density [log(S*) of -1.12 Mq yr"' kpc'^]. 

In terms of environment, we find that 90%-95% of objects 
in the Main, Offset-SF, and Offset-AGN samples have zero 
or one neighbor, the lowest density bin in Kauffmann et al. 
(2004). Furthermore, the distributions of environments for 
the Offset samples are similar to that of the Main sample. 
Therefore, it appears that the vast majority of galaxies consid- 
ered here reside in the lowest density environment defined by 
Kauffmann et al. (2004). This result is not surprising, as our 
SDSS samples all contain emission-line galaxies, which are 
more likely to be found in low-density environments (Kauff- 
mann et al. 2004). Currently we have too few Offset objects 
to draw any solid conclusion about the question on galaxy en- 
vironmental dependence. 

Groves et al. (2006) suggest that the offset on the BPT di- 
agram seen in some high-redshift galaxies is mostly caused 
by contribution from AGNs, as they found properties of their 
candidate low-metallicity H II/AGN composites similar to 
those of the host galaxies of AGNs (Kauffmann et al. 2003a; 
Heckman et al. 2004); these objects have on average signif- 
icantly higher stellar masses, older stellar populations, and 
redder colors than the sample of pure star-forming galax- 
ies used for comparison. Groves et al. (2006) arrive at this 
conclusion after constructing their offset sample in a region 
bounded by the Ka03 and the KeOl curves in the [N II]/Ha 
versus |0 IIIJ/H/j diagram, and the Seyfert branch [log([0 
III]/H/3) > 31og([N II]/Hq)]. When we reproduced their off- 
set sample, we found the median values of log([N II]/Ha) ~ 
-0.286 and log([0 IIIJ/H/3) ~ -0.024, whereas the correspond- 
ing values for our DEEP2 objects are -0.638 and 0.344. As 
the density of objects increases rapidly towards increasing [N 
II]/Hq, the offset sample definition of Groves et al. (2006) is 
weighted heavily towards H II/AGN composites and not nec- 
essarily representative of the average [N IIJ/Ha or [O in]/H/3 
of high-redshift galaxies. Furthermore, it is not clear that their 
"typical" comparison sample, which was defined as all galax- 
ies within ±0.05 dex of |-0.55, 0.10] in the [N II]/Ha versus 
[O III]/H/3 diagram, represents a fair one, as they have not 
controlled for any galaxy property. Given the strong correla- 
tions among galaxy properties, it is crucial to compare galax- 
ies that are similar in at least some basic parameters. Here 
we want to emphasize that in order to make a controlled com- 
parison, the Main sample was selected to have similar range 
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and median of stellar masses as that of the Offset sample. The 
Main sample also has a similar median [N ll]/Ha value. With 
this controlled comparison of typical and Offset SDSS ob- 
jects, we find that the unusual objects on the BPT diagram 
are populated not only by likely AGN hosts (Offset- AGN) but 
also by the Offset-SF sample, most of which do not resemble 
typical AGN-host galaxies. Furthermore, there is no segre- 
gation of the two offset samples in the [N II]/Ha versus [O 
III]/H/3 diagram. 

As suggested in Groves et al. (2006), another possible test 
of the presence of an AGN would require the detection of ei- 
ther He II A4686 or [Ne V] A3426 Hnes. The [Ne V] A3426 
line is not redshifted into the SDSS band for the majority of 
our objects: we have only four such objects in Offset- AGN, of 
which the spectra near [Ne V] A3426 are too noisy to provide 
any meaningful constraints. The He II A4686 line is expected 
to be very weak: a 20% AGN contribution to H/3 implies He 
II/H/3 = 0.05. Note that besides the AGN photoionization, 
possible mechanisms for producing He II emission also in- 
clude hot stellar ionizing continua and shock excitation (Gar- 
nett et al. 1991). Another important goal is to infer metallic- 
ities of the Offset-SF objects using a method that is indepen- 
dent of strong-line indicators, since these objects have abnor- 
mal strong-line ratios. A more calibration-independent way 
to address metallicity is through the classic method, which 
relies on measuring weak auroral lines. Both the test of the 
presence of AGNs, and the determination of metallicity with 
the direct method, require the measurement of very weak 
lines. These measurements become feasible through the use 
of composite spectra. 

5.5. SDSS Composite Spectra 

In order to measure He II A4686 as another test of possible 
AGN contribution to the ionizing spectrum, and [O III] A4363 
to infer the oxygen abundance using the method, we con- 
struct composite spectra for both Offset-SF and Offset- AGN 
objects. For comparison, we also make composite spectra for 
the Main sample. We adopt the same method of making com- 
posite spectra used for our DEEP2 objects, employing median 
scaling to preserve the relative fluxes of the emission features. 
The variance spectrum is combined in the same way and the 



error for the corresponding composite spectrum is the square 
root of the composite variance spectrum divided by the num- 
ber of objects. 

The stellar continuum was modelled over the rest-frame 
wavelength range of 3800 - 7500 and then subtracted from the 
composite spectra, with a method similar to the one described 
in Tremonti et al. (2004) and Brinchmann et al. (2004). This 
technique (Chariot & Bruzual 2007, private communication) 
relies on the Bruzual & Chariot (2003) stellar population syn- 
thesis model, but with two major improvements. First, the 
number of the metallicity grids was expanded from three to 
four; second, and most importantly, the stellar continuum 
fitting was carried out using both the STELIB (Le Borgne 
et al. 2003) and the MILES (Sanchez-Blazquez et al. 2006) 
libraries of observed stellar spectra. Models using the MILES 
library fit the data better than those with the STELIB one, in 
terms of the agreement around Balmer lines. It is especially 
essential to obtain a reasonable stellar continuum fit near H7, 
which is very close in wavelength to the weak [O III] A4363 
auroral line. Models using the STELIB library tend to over- 
estimate stellar continua, especially around Balmer lines. We 
therefore adopt the stellar continuum models based on the 
MILES library. 

The resultant continuum-subtracted spectra, along with the 
original composite spectra and the best-fit model of the stellar 
continua, are shown in Figure (12). We only present portions 
of the entire composite spectra to enable scrutiny of the rel- 
evant weak lines ([O III] A4363 and He II A4686) and our 
capability of making a reasonably good stellar continuum fit. 
The corresponding error spectra are estimated by combining 
the measurement error for the composite spectra and the sys- 
tematic uncertainties from the stellar continuum fitting. The 
systematic error from the stellar continuum fitting is estimated 
using the rms of a portion of the continuum-subtracted spec- 
trum, which is relatively free of emission-line features. This 
method of estimating continuum fit error is motivated by the 
fact that the rms would be zero if the fitting was ideally good 
and if there were no emission-line features. At all of the 
wavelengths considered, the systematic error is significantly 
larger than the measurement error for the composite spectra 
and therefore dominates the total error budget. Emission-line 
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Fig. 12. — SDSS composite spectra. For clarity, only the portion of the spectrum containing the weak emission lines of interest is presented here. The original 
composite and stellar-continuum-subtracted spectra for samples (A) Main, (B) Offset-AGN, and (C) Offset-SF are shown as black solid curves. Stellar-continuum 
fits using Brazual & Chariot (2003) population synthesis models are plotted as red dotted curves. The [O III] A4363 and He II A4686 lines are marked by dotted 
lines. 

fluxes and flux ratios with the total uncertainties are given in 
Table 5. 

There are a few caveats that must be mentioned in the anal- 
ysis of composite spectra described here. First, we are aver- 
aging over tens of thousands of galaxies for the SDSS Main 
sample and dozens of galaxies for the Offset-SF and Offset- 
AGN samples; second, for each single galaxy, we are dealing 
with integrated spectra containing contributions from multiple 
H II regions, which may show a range of metallicities, ion- 
ization parameters, ionizing spectra, and electron densities. 
Also, even though we select the samples to have fiber aper- 
tures covering more than 20% of the total g-band photons, 
because of the presence of radial gradients in galaxy proper- 
ties, uncertainties still remain since the spectrum is weighted 



towards the nucleus. Despite these caveats, however, we want 
to emphasize that the analysis of composite integrated fibre- 
based spectra is still meaningful in terms of determining av- 
erage properties and the relative differences among samples. 

5.5.1. AGN Contribution 

The composite spectra we construct can be used to address 
the level at which AGN ionization contributes to the anoma- 
lous emission-line ratios, discovered in both nearby SDSS 
Offset objects and high-redshift DEEP2 galaxies. Based on 
their [S ll]/Ra and [O l]/Ha ratios, we expect the Offset-AGN 
objects to have a non-negligible contribution to their ioniza- 
tion from AGNs and/or shock heating. Based on the compos- 
ite spectrum, we can also probe the weaker He II line, which 
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TABLE 5 

Emission-line FLUX measurements of SDSS composite spectra. 

Sample F[0IIJAA3726.3729'' F[0IIIJA4363'' FHeIIA4686'' Fh/3'' F[0III]A4£)5£)'' FlOIlIJASOOV'' F[OI]A6300'' 

Main 327.3±0.3 0.6±0.2 1.4±0.2 153.5±0.2 45.2±0.1 137.0±0.3 19.1±0.2 

Offset-AGN 271.9±0.8 2.8±0.5 4.9±0.5 90.4±0.5 86.4±0.2 261.7±0.6 24.3±0.6 

Offset-SF 629.1±1.0 7.1ifc0.6 12.5±0.7 321.0il=0.7 266.6±0.3 807.4±0.8 38.8±0.7 

Sample Fh„- F[sniA67i7- Ftsn]A673i- ^'^S^' '^If^ "[s^amviv^tb. 

Main 633.6±0.3 117.6±0.3 83.9±0.3 249±82 0.009±0.001 0.0301±0.0004 0.3180±0.0006 

Offset-AGN 337.3±0.7 80.6±0.6 58.5±0.6 110±19 0.054±0.006 0.0720±0.0018 0.4124±0.0028 

Offset-SF 1272.8±0.9 170.6±0.8 137.3±0.8 129±11 0.039±0.002 0.0305±0.0006 0.2419±0.0009 



^Emission-line flux and 1 a error in units of 10 ergs s ^ cm ^ from SDSS composite spectra. 
^Dust-reddening corrected using Balmer decrement method. 



is another indicator of AGN activity. 

From the He II/H/3 ratio (0.054±0.006), we can see that 
the Offset-AGN sample may have up to ~20% AGN con- 
tribution to the Balmer emission lines (Groves et al. 2006). 
The Offset-SF sample, on the other hand, shows a lower level 
of possible AGN contamination (He ll/H/3 of 0.039±0.002). 
This distinction is also suggested by the different host-galaxy 
properties of the two samples. The Offset-AGN objects with 
similar stellar masses to the Main sample have redder colors, 
older stellar populations, lower starburst stellar mass fraction, 
smaller Ha EWs, and smaller specific star formation rates, 
which are similar to host galaxies of AGN (Kauffmann et al. 
2003a), whereas the Offset-SF sample is similar to the Main 
sample in terms of these galaxy properties. Furthermore, by 
construction, the Offset-SF objects have low [S II]/Ha and [O 
1]/Hq; in the range favored for typical star-forming galaxies. 

We note, however, that the Offset-SF sample still 
has a higher value of He II/H/3 than the Main sample 
(0.009±0.001). There might be some mixing with the Offset- 
AGN objects because of uncertainties in emission-line diag- 
nostic ratios. Indeed, three objects in the Offset-SF sam- 
ple have cross-identifications with ROSAT sources, and two 
of them have optical spectra indicating broad Balmer fines. 
However, ~90% of the Offset-SF sample show no evidence 
of AGNs in terms of ROSAT-somce cross-identifications and 
broad Balmer lines. In addition, we have examined the He 
Il/H/3 ratios from individual galaxy spectra in the Offset-SF 
sample, finding that five objects (including the three with 
ROSAT cross-identifications) of the 25 have He II/H/3 ~ 0.05- 
0.10, while the remaining 20 objects all have He II/H/3 < 0.01, 
similar to the typical He II/H/3 of the Main sample. Even 
the five objects with He Il/H/3 0.05-0.10 could have other 
possible ionizing sources including hot stars and shocks (Gar- 
nett et al. 1991). It is worth mentioning that nondetection in 
ROSAT or broad Balmer fines cannot fully exclude the pos- 
sibility of AGN contamination, as the majority of obscured 
AGNs are completely absorbed in the soft X-ray range and 
do not show broad Balmer fines. However, if the dominant 
cause of the Offset-SF on the BPT diagram for our Offset- 
SF sample was mildly obscured AGN contamination, namely, 
if their optical spectra were contaminated by AGN ioniza- 
tion at a lower level than the Offset-AGN sample, then we 
would expect them to exhibit intermediate emission-line diag- 
nostic ratios and host-galaxy properties, relative to the Main 
and the Offset-AGN samples. This is not the case. Instead, 
the Offset-SF objects have even higher ionization parameters 
and electron densities than the Offset-AGN sample, relative 
to the Main sample. In addition, as discussed in §5.4, Offset- 



SF objects have host-galaxy properties different from those 
of typical weak AGNs in SDSS (Kauffmann et al. 2003a). 
Therefore, it appears that AGN excitation does not provide 
the dominant cause of the anomalous emission-fine ratios for 
most of the Offset-SF sample. 

5.5.2. Electron Temperature and Oxygen Abundance 

By making composite spectra, and hence enhancing the S/N 
of weak auroral lines, we can determine metallicities with a 
method independent of strong-line indicators. We have al- 
ready discussed the fact that the Offset-SF objects have larger 
ionization parameters, as indicated by their higher val- 
ues. However, it is still unclear whether the Offset-SF sample 
has comparable metallicity with the Main one, on average, as 
suggested by N2 and N202, or if it has ^0.3 dex lower metal- 
licity, as suggested by 03N2 and 7?23- Settling this question is 
important for multiple reasons. 

First, gas-phase metallicity and electron temperature, them- 
selves, serve as key parameters of H II regions. Second, 
understanding the relative metallicities of the Main and the 
Offset-SF samples on average can help us determine the dif- 
ference in average ionization parameter quantitatively. If the 
Offset-SF sample has comparable metalficities with the Main 
sample, their ionization-parameter difference would be very 
large; on the other hand, if the Offset-SF sample has much 
lower metallicities than the Main one does, their ionization- 
parameter difference would be much smaller, since both lower 
metalficity and higher ionization parameter can cause large 
O32. Finally, if the Offset-SF objects actually have much 
lower metallicities than the Main objects, then they will not 
follow the local mass-metallicity relation (Tremonti et al. 
2004), given their comparable stellar masses. 

In order to determine oxygen abundances using the direct 
method, estimates of both the electron density and temper- 
ature are required. Based on the ratio [S II] A6717/[S II] 
A6731, we establish that on average all of our objects are 
in the low-density regime. Then using the five-level atom 
program nebular implemented in IRAF/STSDAS (Shaw & 
Dufour 1995), we derive the electron temperature T[0 III] 
from the ratio of auroral line [O III] A4363 to nebular lines [O 
111] AA4959, 5007 (Osterbrock 1989) and correct for dust ex- 
tinction using the Balmer decrement method. As for T{0 II], 
we adopt a simple scaling relation between the temperatures 
in different ionization zones of an H 11 region predicted by the 
photoionization models of Garnett (1992): 

r[011]=0.70r[01II] + 3000K, (3) 

which is applicable in a wide range of r[0 III] (2000-18,000 
K). loruc abundances O^IR^ and 0"^/H"^ are then determined 
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Fig. 13. — Mass-metallicity relation for the SDSS Main (grey contours) and Offset-SF (red open circles) samples, based on various strong-line indicators. 
Calibrations are from Zaritsky et al. (1994, for R23), Pettini & Pagel (2004, for N2 and 03N2), and Bresolin (2007, for N202), respectively. 



TABLE 6 

Physical Quantities from SDSS composite spectra. 



Sample N/ T^LOIII]'' Tf[On]'' 0+^+"^^ 0++^+"^^ 12 + log(0/U)'^ 

Main 28±5 0.950+0 '2;* '^■'^^^-i'ofo ^^■^'^'] 3.7±1.3 8.23+j; |l 

Offset-SF 190±12 1-170™^ 1119!!J;|]?;^ 6.5±0:6 5.2±0.5 8.07±0.04 



^Electron density and 1 cr error in units of cm ^ derived from ([S II]A6717/[S I1JA673I). 
^Electron temperature and 1 cr eiTor in units of 10'* K. 

'^lonic abundance and 1 a eiTor from coUisionally excited lines in units of 10"^. 
'^Total oxygen abundance and I a enor assuming O/H = (0+ + 0^)/H+. 



using programs in the nebular package, and, finally, we 
obtain the total gas phase oxygen abundance assuming O/H = 
(O""" + 0^)flV'. The results are summarized in Table 6. 

First, the Offset-SF sample has ^2200 K higher electron 
temperatures than the Main one does on average. These ob- 
jects are unusual in the sense that they have higher electron 
densities and electron temperatures, hence ambient interstel- 
lar pressures assuming pressure equilibrium between H II re- 
gions and ambient gas, and larger ionization parameters. It 
is not clear what is at the root of the special H II region 
physical conditions, in terms of galaxy properties, or larger 
scale environments. However, in the simple theoretical es- 
timates for local starburst galaxies, the ISM pressure scales 
roughly linearly with the SFR surface density (e.g., Thomp- 
son et al. 2005). Indeed, the Offset-SF sample has several 
prominently related properties, including higher concentra- 
tion (r-band Petrosian concentration index 7?9o//?5() of ~2.6 
compared to ~2.4), smaller half-light radii [.R5o(z) of 1.2 kpc 
compared to 1.9 kpc], and most notably, higher SFR surface 
density [log(S*) of -0.86 M0 yr"' kpc"- compared to -1.12 
Mq yr"' kpc"^)], as shown in Figure 10, compared to typical 
star-forming galaxies. This higher SFR surface density may 
therefore account for the higher interstellar pressure seen in 
the H II regions of Offset-SF objects. 

It is also possible to relate electron temperature, electron 
density, and ionization parameter in a single H II region, 
through the classical picture of Stromgren (1939). With 
higher ambient interstellar pressure due to higher electron 
densities and temperatures, H II regions are surrounded by 
denser molecular gas dust and therefore have smaller radii un- 
der the stall condition. In an idealized case of a fully ionized 
spherical H II region with pure hydrogen, the radius R of the 
ionized region is given by R= [3Q/(4-nN^aB)] where Q is 
the rate of emission of hydrogen ionizing photons, and ag is 
the case B recombination coefficient (Osterbrock 1989). The 
ionization parameter U defined by [F /AttcR^Nhu], where F 



is the flux of ionizing photons, c is the speed of light, and A^hii 
is the particle density in H II region, is then proportional to 

1/3 

A^;'-'for a given central hot star Therefore, in this simple pic- 
ture, H II regions with higher electron densities and temper- 
atures will have higher ionization parameters. Accordingly, 
the higher electron densities and temperatures inferred from 
the integrated spectra of the Offset-SF objects may result in 
the observed higher ionization parameters as well. 

Second, the Offset-SF sample has ^0.16 dex lower metal- 
licities than the Main sample. Note that the absolute abun- 
dance values inferred using the method are significantly 
lower than those based on any strong-line indicators. This 
difference may stem from temperature fluctuations in H II re- 
gions and the strong Tg dependence of line emissivities, from 
which the net effect is that the electron temperature derived 
from auroral lines tends to overestimate the real temperature 
and hence systematically underestimate the abundance (Gar- 
nett 1992). Therefore, we only focus on relative differences 
based on the same method between the Offset-SF and Main 
samples. Given that the Offset-SF and Main samples are char- 
acterized by similar average stellar masses, this difference 
in average metallicity suggests that Offset-SF objects do not 
follow the local mass-metallicity relation, in addition to be- 
ing more compact and sustaining higher interstellar pressures 
than the Main sample. 

To investigate this question further, we examine several ver- 
sions of the mass-metallicity relation for the Offset-SF and the 
Main objects, based on different strong-line metallicity indi- 
cators. Figure 13 includes these mass-metallicity correlations 
for both the Main and the Offset-SF samples, demonstrating 
that one would draw different conclusions about the relation- 
ship between the Offset-SF and the Main samples, depending 
on which strong-line indicator was used. Between the Main 
and Offset-SF samples, the relative average metallicity differ- 
ences at fixed mass based on 7?23, N2, 03N2, and N202 are 
^^0.28, 0.00, 0.18, and 0.00 dex, respectively. This ambiguity 



20 



LIU ET AL. 



could be understood as the result of adopting the same calibra- 
tion for samples with different physical conditions, of which 
the most important are ionization parameter and electron den- 
sity. 

The discrepancy between conclusions drawn from the 
strong-Une indicators and the Te method can be understood 
using the photoionization models of Kewley & Dopita (2002). 
For example, according to Rxi, the average metallicity differ- 
ence between the Offset-SF and Main samples is ^0.28 dex, 
whereas this metallicity difference based on the Tg method 
is ~0.16 dex. This ~0.12 dex discrepancy can be explained 
by examining the theoretical grids in the Rxi versus 12 H- 
log(0/H) space, given by photoionization models with differ- 
ent values of electron density and ionization parameter. If the 
incorrect assumption is adopted that both the Main and Offset- 
SF samples have electron densities A'e ~ 10 cm"' and ioniza- 
tion parameters = 3 x lO^cm s"' - values appropriate for the 
Main but not Offset-SF sample - we obtain that the Offset- 
SF sample is ~ 0.25 dex lower in oxygen abundance than 
the Main sample. On the other hand, if an electron density 
Ne ~ 350 cm~^ and ionization parameter ^ = 5 x lO^cm s~^ are 
assumed for the Offset-SF sample - more appropriate based 
on their inferred H II region physical conditions - the differ- 
ence in estimated Main and Offset-SF oxygen abundance is 
only ~0.10 dex. This difference is ~0.15 dex smaller than 
the one inferred with the incorrect assumption that the Offset- 
SF and Main samples have the same electron densities and 
ionization parameters. The actual difference in electron den- 
sity between the two samples is less extreme (A^^ ~ 190 cm"' 
of the Offset-SF compared to Ne ~ 28 cm"' of the Main), so 
this discrepancy may be smaller, as observed (~0.12 dex). 
Analogous effects occur for the N2 and 03N2 indicators. 

In summary, the metallicities of the Offset-SF objects based 
on strong-line indicators can be either overestimated or under- 
estimated, depending on the relative shifts of the theoretical 
grids caused by varying physical conditions, and an incorrect 
assumption of the values of ionization parameter and electron 
density. For the SDSS Offset-SF sample, the metallicities 
based on strong-line indicators can either be systematically 
too large by ~0.16 dex (N2, N202), or too low by ~0.12 dex 
(/?23). The bias of 03N2 is little (too large by ~0.02 dex) for 
our SDSS Offset-SF objects, although it may not always be 
negligible given other ranges of physical conditions in terms 
of metallicity, ionization parameter, and electron density. The 
reason for the Offset-SF objects' anomalous positions on the 
BPT diagram is the same: H II regions with distinct ionization 
parameters and electron densities fall on different surfaces in 
the diagnostic-line parameter space. 

One caveat that should be mentioned here is the fact that we 
are inferring average biases over the stellar mass and metalUc- 
ity ranges spanned by the Offset-SF sample. Another caveat is 
that results of the SDSS data are all based on integrated fiber 
spectra and hence limited in terms of spatial information. We 
plan to further study the spatial dependence of H II region 
emission Unes in these offset objects using long-sUt spectro- 
graphs in the future. 

5.6. Implications for DEEP2 Objects 

We have analyzed emission-line diagnostic ratios, physical 
conditions and galaxy properties of SDSS objects with similar 
[N II]/Ha and [O III]/H/3 values to those of DEEP2 galaxies 
with the most extreme offset in our sample, and z ~ 2 star- 
forming galaxies (Erb et al. 2006a). There are two major 
causes for their offset: one is different H II region physical 



conditions characterized by higher electron density and tem- 
perature, and hence larger ionization parameter, compared to 
typical SDSS star-forming galaxies. These physical condi- 
tions of the SDSS Offset-SF objects are also coimected to 
their host galaxy properties, particularly the higher SFR sur- 
face density. The other possible cause for the offset is con- 
tribution from AGNs and/or shock excitation. We cannot rule 
out either of these two possibilities, or a combination of both, 
to explain the offset of DEEP2 objects on the diagnostic dia- 
gram. 

As for the question of H II region physical conditions, it 
is possible to determine electron densities for our DEEP2 ob- 
jects using the [O IIJ doublet contained in the DEIMOS spec- 
tra. Note that we will not directly compare electron densities 
of the z ~ 1.0-1.5 objects to those of SDSS local galaxies, 
because they were estimated using different density-sensitive 
doublets, and therefore the comparison might suffer from the 
associated systematic uncertainties. However, when we di- 
vide all of the objects in our DEEP2 sample that have reliable 
[O II] -doublet measurements into two groups, according to 
the empirical Ka03 curve in the [N II]/Ha versus [O III]/H/3 
diagram, we find that for the group below the Kauffmann 
curve, the median value of electron density inferred from [O 
II] A3726/[0 II] A3729 is ~23 cm"', whereas for the group 
above the Kauffmann curve, the median value is ~159 cm"'. 
This difference in electron density as a function of position 
on the BPT diagram is qualitatively consistent with what we 
have found for SDSS objects (Table (6); electron density Ng 
[cm"'] of 190±12 for the Offset-SF compared to 28±5 for 
the Main), lending independent support for our method of us- 
ing Offset-SF objects as local analogs for the H II regions in 
our high-redshift sample. In principle, we could examine the 
ionization parameter indicator, O32, for our DEEP2 objects, 
using [O H] from DEIMOS and [O IH] from NIRSPEC spec- 
tra. However, at this point, the comparison does not seem 
useful, as it would suffer from large systematic uncertainties, 
due to the manner in which the optical and near-IR spectra 
were collected and flux-calibrated. 

A related question is whether objects that are more offset in 
the [N II]/Ha versus [O III]/H/? diagram have higher average 
SFR surface density than objects with less offset within our 
DEEP2 z ^ 1-1.5 sample. The spatial extent of Ha emis- 
sion contains information about galaxy sizes. However, we 
can only measure the extent along the slit, which does not 
necessarily reflect the actual galaxy size. In addition, the size 
estimated from the spatial extent of Ha emission is subject to 
uncertainties due to the variation in the seeing FWHM, which 
is comparable to the Ha size itself. As we do not have reli- 
able estimates for the sizes of our DEEP2 objects, the above 
question can be finally answered only when the galaxy mor- 
phological information is robustly gathered. Furthermore, the 
SFRs of our DEEP2 objects have not been corrected for dust 
extinction or aperture effects, which may amount to factor of 
2 differences (Erb et al. 2006b). Limited by the associated 
systematic uncertainties, we are not able to compare the SFR 
surface densities of our DEEP2 objects directly with those 
of the SDSS samples. However, after correcting the SFR- 
associated systematic uncertainties for their z ~ 2 sample, 
Erb et al. (2006c) found a mean log(S,t) « 0.46, significantly 
higher than that of SDSS local star-forming galaxies. In ad- 
dition, observational evidence exists that the galaxy size at 
fixed mass/luminosity decreases with increasing redshift out 
to z ~ 3 (Trujillo et al. 2006; Dahlen et al. 2007). At the same 
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time, galaxies with star formation rates significantly higher 
than those found among typical SDSS emission-line galaxies 
are more commonly found at higher redshifts (e.g., Dahlen 
at al. 2007; Reddy at al. 2007). Tharafora, galaxies with high 
SFR surface densities are more prevalent in the high-redshift 
universe. For such objects, a noticeable difference in H II re- 
gion physical conditions is expected. 

As for the question of AGN contribution, the DEEP2 ob- 
ject 42010637 clearly falls on the part where many objects 
have contribution from AGNs in their ionizing spectra. How- 
ever, it does not have multi-wavelength information that can 
either confirm or rule out the AGN excitation. In addition, 
all of the z ~ 2 objects with [O m]/H/3 and [N II]/Ha mea- 
surements are offset by an amount similar to our most ex- 
tremely offset z ~ 1.0-1.5 objects but show no evidence of 
AGN contamination in at least their rest-frame UV spectra 
(Erb et al. 2006a). The lack of such evidence in the rest- 
frame UV rules out the presence of at least a fairly unobscured 
AGN. Note that Daddi et al. (2007) find that roughly 20-30% 
of star-forming galaxies with M ~ lO'^-lO'^M© at z ~ 2 dis- 
play a mid-IR excess, as evidence for hosting an obscured 
AGN. These authors show that this fraction increases with 
stellar mass, reaching ^^50-60% for an extreme mass range 
of M > 4 X 1O"'M0. Our DEEP2 z ~ 1.0- 1.5 sample covers 
a range of M ~ 5 x 10^- 10' 'M©, not all of which are in the 
"massive" range of Daddi et al. (2007), and only two of 20 
are in the extreme range. If the mid-IR excess, hence AGN 
contamination, found by Daddi et al. (2007) is also applicable 
to our DEEP2 z ~ 1.0-1.5 objects, the fraction of objects 
containing potential AGN contamination would be at most 
~ 20%-30%. However, nine (or 10, one with a lower limit) of 
the 13 objects in our sample with all four line measurements 
show evidence for being offset in the [O III]/H/3 versus [N 

II] /Ha BPT diagram (i.e., more than 20-30%). More gener- 
ally, the associated space densities of both the blue DEEP2 
galaxies (Coil et al. 2008) and UV-selected z ~ 2 objects for 
which we have NIRSPEC spectra (Adelberger et al. 2005; 
Reddy et al. 2007) are significantly higher than that of the 
obscured AGN population featured in Daddi et al. (2007). 
Therefore, as one of the possible causes for the offset on the 
BPT diagram, AGN contamination cannot account for nor is 
consistent with all the rest-frame optical emission-line mea- 
surements presented here. 

We will address the AGN-contamination issue in the fu- 
ture by looking at more DEEP2 objects for which multi- 
wavelength information is available and obtaining spatially 
resolved spectra with an integral field unit assisted by adap- 
tive optics to isolate the contribution from the nucleus. Also, 
the measurements of flux-calibrated emission lines including 
[O II] A3727, [O III] A4363, H/3, [O III] A5007, [O I] A6300, 
[N n] A6584, Ha and [S II] AA6717, 6731, as well as host- 
galaxy morphological information are required, in order to fi- 
nally settle the causes of the offset in the [N II]/Ha versus [O 

III] /H/3 diagram for high-redshift star-forming galaxies. [O 
III] A4363 may be difficult to measure for all but the most 
metal-poor objects or in deep z ~ 1 composite spectra, and 
intensity limits on [O I] A6300 may also be difficult to ob- 
tain. Yet measurements of the stronger emission lines for a 
statistical sample of objects at z > 1 will be feasible with the 
next generation of ground-based multi object near-IR spectro- 
graphs. 

In an independent study, Brinchmann et al. (2008) have also 
analyzed the possible causes for the high-redshift galaxies' 
offset in the BPT diagram, using theoretical models of nebu- 



lar emission from star-forming galaxies (Chariot & Longhetti 
2001) and the SDSS DR4 data. They have found a relation- 
ship in SDSS galaxies between their location in the BPT di- 
agram and their excess specific SFRs and larger Ha EWs 
relative to galaxies of similar mass. We note that they have 
only examined SDSS star-forming galaxies below the Ka03 
curve, whereas our SDSS Offset samples have been selected 
to be above this curve, in order to probe exactly the regime 
where the most offset objects in our DEEP2 z ~ 1 - 1 .5 sam- 
ple and the Erb et al. (2006a) z-2 star-forming galaxies 
reside. Brinchmann et al. (2008) have inferred that an el- 
evated ionization parameter U is at the root of the excess 
specific SFRs of the more offset objects within their pure 
star-forming galaxy sample and further speculated that higher 
electron densities and escape fractions of hydrogen ionizing 
photons might be the factors responsible for the systemati- 
cally higher values of U in the H II regions of high-redshift 
galaxies. Using a different technique and sample of galaxies, 
we have reached a similar conclusion about the higher ion- 
ization parameter and larger Ha EWs in our SDSS Offset- 
SF samples. In addition, we have also uncovered that the 
higher electron density and temperature, hence higher inter- 
stellar ambient pressure, is at the root of the higher ionization 
parameter. We have further shown that these unusual H II re- 
gion physical conditions are well connected to the higher SFR 
surface density of host galaxies. The trend of higher electron 
density with increasing BPT diagram offset found within our 
DEEP2 z 1 - 1 .5 sample, and the observational evidence that 
galaxies with high SFR surface densities are more prevalent at 
high redshifts, lend further support to our conclusions drawn 
based on the SDSS local emission-line galaxies. 

These differences discovered in H II region physical condi- 
tions, which may commonly apply to z ^ 1-1.5 star-forming 
galaxies, must be taken into account when strong-line abun- 
dance indicators are used to study the evolution of galaxy 
metallicity with redshift. The resulting systematic bias in 
inferred oxygen abundance can be estimated quantitatively 
either via detailed photoionization models, given the differ- 
ence in H II region physical conditions inferred from the rele- 
vant density- and ionization parameter-sensitive line ratios, or 
through empirical comparisons, as we have illustrated for the 
SDSS Main and Offset-SF samples. 

We have shown in §5.5.2 that, for the SDSS Offset-SF ob- 
jects, the metallicities based on various strong-line indica- 
tors can either be systematically too large by ~0.16 dex (N2, 
N202) or too low by ~0.12 dex (/?23). For N2 in particular, 
the inferred metallicities of the SDSS Offset-SF objects can 
be systematically too large by ^0.16 dex, which is already 
comparable to the inherent scatter in the N2 calibration ( 1 a 
dispersion of 0.18 dex; Pettini & Pagel 2004). This system- 
atic bias may be even larger if the actual difference in H II 
region physical conditions is more extreme. In addition, we 
note that this systematic uncertainty stemming from the dif- 
ference in H 11 region physical conditions has different effects 
from that of the inherent scatter in the caUbration, when the 
strong-line relation calibrated with local H II regions is ap- 
pHed to high-redshift star-forming galaxies. As there is evi- 
dence that the average ionization parameter and electron den- 
sity in high-redshift star-forming galaxies are systematically 
higher than the local typical values, there will be a system- 
atic "bias," instead of a scatter as a result. For example, the 
N2-based metallicities for the DEEP2 z ^ 1.4 sample may be 
systematically overestimated by as much as ^ 0.16 dex. Fur- 
thermore, as discussed in §4, the fact, that our DEEP2 z ^ 1.4 
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sample is more offset from the local excitation sequence than 
the z '--^ 1 .0 sample, may result in N2-based metallicities that 
are more significantly overestimated as well. This effect may 
cause the apparent reverse trend of average O/H with redshift 
within the DEEP2 sample, such that the z ~ 1 .4 sample is de- 
scribed by apparently higher metallicities at fixed stellar mass 
than the one at z ^ 1.0, despite the general trend of increas- 
ing O/H towards lower redshift. The reverse trend within the 
DEEP2 sample based on 03N2 is less significant than that 
based on N2, which is consistent with the fact that the 03N2- 
based metallicity bias due to the offset in the BPT diagram is 
much less than that of N2. 

In addition, we can quantify the potential bias for the sam- 
ple of star-forming galaxies at z ~ 2 presented by Erb et al. 
(2006a). If the four z - 2 objects with the full set of H/3, [O 
III] A5007, Hq, and |N II] A6584 measurements are represen- 
tative of the larger population of star-forming galaxies in Erb 
et al. (2006a) in terms of diagnostic line ratios, such that the 
z ~ 2 star-forming galaxies are even more significantly offset 
than the DEEP2 z 1.4 sample on average, then their metal- 
licities based on N2 would also be overestimated by ~0.16 
dex. Therefore, the true metallicity offset between the Erb 
et al. (2006a) z ~ 2 sample and local SDSS objects may be 
as much as ~0.16 dex larger than it appears now in Figure 
(6a). Accounting for the systematic differences in converting 
strong emission-line ratios to oxygen abundances is therefore 
a crucial component of comparing galaxy metallicities at dif- 
ferent redshifts. 

6. SUMMARY 

We have compiled a sample of 20 star-forming galaxies at 
1 .0 < z < 1.5 drawn from the blue cloud of the color bimodal- 
ity observed in the DEEP2 survey, to study the correlation be- 
tween stellar mass and metallicity, across a dynamical range 
of 2 orders of magnitude in stellar mass, as well as H II re- 
gion physical conditions at this redshift range. In order to gain 
some insights on the causes of the offset in the BPT diagram 
observed in high-redshift star-forming galaxies, we have ex- 
amined the H II region diagnostic line ratios and host galaxy 
properties of the small fraction of SDSS galaxies that have 
similar diagnostic ratios to those of the DEEP2 sample. Our 
main results are summarized as follows: 

1. There is a correlation between stellar mass and gas- 
phase oxygen abundance in DEEP2 star-forming galax- 
ies at z ^ 1.0 and at z ^ 1 .4. We have found that the 
zero point of the M-Z relationship evolves with red- 
shift, in the sense that galaxies at fixed stellar mass be- 
come more metal-rich at lower redshift, by comparing 
the 1.0 < z < 1.5 sample with UV-selected z ~ 2 and 
SDSS local star-forming galaxies. At the low-mass end 
(M,^ ~ 8 X IO'^Mq), the relation at 1.0 < z < 1.5 is off- 
set by ~0.2 (0.35) dex from the local mass-metallicity 
relation according to the N2 (03N2) indicator. The N2- 
based offset could be larger by as much as ~0.16 dex, 
when the systematic bias due to difference in H II re- 
gion physical conditions between 1 .0 < z < 1.5 and the 
local universe is taken into account. At the high-mass 
end {Mi, ~ 5 x IO^^Mq), the metalUcity offset between 
the DEEP2 1.0 < z < 1.5 sample and the local SDSS 
sample is at least ~ 0.2 dex, according to the 03N2 in- 
dicator. 

2. As observed previously for a very small sample of high- 
redshift galaxies, on average our new DEEP2 sample 



at 1 .0 < z < 1 .5 is offset from the excitation sequence 
formed by nearby H II regions and SDSS emission- 
line galaxies. By examining the small fraction of SDSS 
galaxies that have similar diagnostic ratios to those of 
the DEEP2 sample, we have found two Hkely causes for 
the anomalous emission-line ratios. One is the contri- 
bution from AGN and/or shock excitation at the level of 
^ 20%. The other is the difference in H II region physi- 
cal conditions, characterized by significantly larger ion- 
ization parameters, as a result of higher electron den- 
sities and temperatures, and hence higher interstellar 
ambient pressure, than the typical values of local star- 
forming galaxies with similar stellar mass. These un- 
usual physical conditions are possibly connected to the 
host-galaxy properties, most importantly smaller sizes 
and higher star-formation rate surface densities. Our 
conclusion drawn from analyzing the SDSS data has 
been further verified by the fact that the DEEP2 objects 
more offset from the local excitation sequence in the 
BPT diagram also have higher electron densities than 
those closer to the local sequence. We cannot rule out 
either the contribution from AGN and/or shock exci- 
tation, or the difference in H II region physical condi- 
tions, for the unusual emission-line diagnostic ratios of 
high-redshift star-forming galaxies. 

3. We have quantified the effects of different H II region 
physical conditions on the strong-line metallicity cali- 
brations. The direct electron temperature method was 
used to estimate the "true" metallicity difference be- 
tween offset SDSS objects with anomalous line ratios 
and more typical objects of similar stellar mass. Strong- 
line indicators were also used to estimate this differ- 
ence. A comparison of these results reveals poten- 
tial biases in the strong-line indicators. According to 
our test, the metallicities based on strong-line indica- 
tors can either be systematically too large by ~0.16 dex 
(N2, N202), or too low by ~0.12 dex (/?23), for objects 
with similar H II region physical conditions to those 
observed in high-redshift galaxies. The bias of 03N2 is 
much less significant (too large by ~0.02 dex) for off- 
set SDSS objects with anomalous line ratios, although 
it may not always be negligible given other ranges of 
physical conditions in terms of metallicity, ionization 
parameter, and electron density. 

The difference in H II region physical conditions, which 
may commonly apply to z '-^ 1 .0 - 1 .5 star-forming galaxies, 
must be taken into account when strong-line abundance in- 
dicators are used to study the evolution of galaxy metallicity 
with redshift. There are at least two methods to remove the 
systematic bias from the effect of significantly different H II 
region physical conditions on the strong-line abundance cal- 
ibrations. One is to gather the abundance information with 
direct Te method for a sample of high-redshift H II regions 
as the calibration sample, which may be hard to achieve, as 
auroral lines are difficult to measure except for very metal- 
poor objects or in deep, composite spectra. The other, which 
is currently feasible yet relies on photoionization models, is 
to quantify the biases in strong-line indicators when certain 
physical conditions are present and then compensate the bi- 
ases when inferring abundances from strong-line indicators 
for high-redshift galaxies. 

In this study we have presented evidence that high-redshift 
star-forming galaxies possess distinct H II region physical 
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properties, as characterized by on average larger ioniza- 
tion parameters, higher electron densities, and temperatures, 
which are possibly connected to their relatively smaller sizes 
and higher SFR surface densities. These conditions may be 
quite common during the epoch at z > 1 when at least 50% of 
the local stellar mass density was formed (Bundy et al. 2006; 
Drory et al. 2005). Therefore, they should be characterized 
in more detail for a full understanding of the star formation 
history of the universe as well as the buildup of heavy ele- 
ments in galaxies. The next generation of ground-based near- 
IR multi-object spectrographs will play a key role in assem- 
bling rest-frame optical emission-line measurements for large 
samples of high-redshift galaxies, enabling the detailed study 
of star-forming galaxies in the early universe. 
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